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ABSTRACT 


A method is presented to integrate the design space of structural/control 
system optimization problems in the case of linear state feedback control. Con- 
ventional structural sizing variables and elements of the feedback gain matrix are 
both treated as strictly independent design variables in optimization by extending 
design variable linking concepts to the control gains. Several approximation 
concepts including new control design variable linking schemes are used to for- 
mulate the integrated structural/control optimization problem as a sequence of 
explicit nonlinear mathematical programming problems. Examples which involve 
a variety of behavior constraints, including constraints on dynamic stability, 
damped frequencies, control effort, peak transient displacement, acceleration and 
control force limits, are effectively solved by using the method presented. 
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Chapter I 

INTRODUCTION 


1.1 INTRODUCTION 

Large space structures usually have low stiffness and low damping char- 
acteristics due to their light weight requirements. In order to suppress the vi- 
bration and maintain the strict shape specifications, it is necessary to enhance 
stiffness and/or damping of the structures through some type of active controls 
(Ref. 1). 

Conventionally structural and controller design is performed separately 
(namely, the structure is designed first to minimize the weight satisfying 
structural constraints and then a control law is found for the fixed structure 
minimizing some kind of control performance index), and as a result the final 
design cannot ensure the best performance of the overall system since the dy- 
namic interaction between the two systems is not directly considered in the 
design process. In Ref. 2 it was shown that slight structural modification can 
lead to considerable improvement in the control system performance, and there 
has been a growing effort to integrate the design optimization of structures and 
control systems in order to achieve a better performance and directly handle 
cross coupling effects and dynamic interactions between the two systems. 


1 


Most of this research has focused on linear control laws based on output 
feedback or state feedback. In the case of output feedback several studies have 
been made where the structural dimensions and the control gains are treated 
as strictly independent design variables in optimization (Refs. 2-10). On the 
other hand, in the case of full state feedback control a sequential approach is 
usually adopted in which the control gains are determined by solving Riccati 
equations corresponding to the changing structural system during design iter- 
ations (Refs. 11-17). When the gain variables are determined by solving 
Riccati equations for a fixed plant, they implicitly become dependent design 
variables and the resulting design optimization is constrained to a subspace 
where the optimality conditions of a control subproblem are satisfied. The 
tendency to subordinate gains to a dependent variable status can be attributed 
to the fact that for system models with a large number of degrees of freedom, 
the gain matrix [//] contains prohibitively large numbers of independent de- 
sign variables (i.e. M x IN control design variables, where M is the number of 
actuators and N is the number of degrees of freedom in the structural model). 

The main purpose of this study is to suggest a new simultaneous opti- 
mization method where both structural and control system variables are 
treated as strictly independent design variables in the case of state feedback 
control. This is accomplished via an adaptation of the design variable linking 
idea to the control system design variables (see Refs. 26 and 27). This inno- 
vation makes it possible to simultaneously optimize the structure/control sys- 
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tem, while avoiding a prohibitively large increase in the total number of 
independent design variables. 


1.2 BACKGROUND 

In Refs. 2-10 both structural and control variables are treated as inde- 
pendent in optimization. Output feedback is adopted as a control law so the 
number of elements in the feedback gain matrix is relatively small. As a result 
the gain elements along with structural design variables can be directly treated 
as independent design variables. 

In Ref. 2 collocated direct velocity output feedback, which is similar to a 
viscous dashpot, is used as a controller. The viscous damping coefficient is 
minimized, with constraints on the closed-loop eigenvalues, by allowing small 
changes of the structural dimensions. 

In Refs. 3-5 in addition to structural dimensions and control gains, 
sensor/actuator locations are also used as design variables. Homotopy and se- 
quential linear programming algorithms are used to optimize either structural 
mass, robustness or eigenvalue sensitivity. Also in Ref. 3 the state feedback 
control case is cast into a similar simultaneous optimization form by using el- 
ements of weighting matrices of the LQR (linear quadratic regulator) as inde- 
pendent design variables, but only the output feedback case is illustrated by 
giving a numerical example. 
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In Refs. 6 and 7 collocated output feedback is chosen as a controller to 
optimize a control augmented structure. In addition to the structural sizing 
variables and control feedback gains, lumped nonstructural masses are treated 
simultaneously as independent design variables. Harmonic dynamic loadings 
are applied and a variety of constraints are considered including natural fre- 
quencies, static displacement and stress, dynamic displacements and actuator 
forces. Due to the characteristics of the collocated sensor/actuator pairs, sys- 
tem stability can be ensured by imposing side constraints on the control gains. 

In Refs. 8 and 9 noncollocated output feedback is chosen instead of 
collocated scnsor/'actuator pairs, and constraints on the stability (real parts of 
closed-loop eigenvalues) are included. 

In Ref. 10 several output feedback control laws are used in the case of 
stochastic disturbances with constraints on the allowable mean square de- 
flection or control effort. 

In Refs. 1 1-25 LQR theory is used for the case of state feedback control. 
In the LQR problem once the weighting matrices in the quadratic performance 
index are chosen, the control law (or all the closed-loop characteristics) is de- 
termined by solving a nonlinear matrix Riccati equation. So the choice of 
weighting matrices in the LQR problem is very important and two problems 
arise. One is how to select meaningful weighting matrices, and the other is the 
solution of the Riccati equation for problems involving a large number of de- 
grees of freedom. 
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In Ref. 12 structural mass is minimized using structural design variables 
while satisfying open-loop frequency constraints and then the LQR problem 
is solved for the fixed structure with given weighting matrices. Here weighting 
matrices are chosen such that the quadratic performance index represents the 
absolute weighted sum of kinetic, strain and potential energies, and the effect 
of relative weighting of these energy terms is discussed. In Refs. 11 and 13 
structural variables are optimized with constraints on the closed-loop 
eigenvalue and modal damping ratios, then the LQR problem is solved for the 
fixed structure with given weighting matrices (identity matrices in this case). 
In Refs. 14-17 the Frobenious norm of the gain matrix is introduced as either 
an objective or a constraint. 

Refence 20 points out the difficulties of simultaneous structural/control 
design and suggests optimization of the closed-loop system using only struc- 
tural tailoring. In this case the objective of structural tailoring is to maximize 
modal stiffness in order to minimize control effort. The control law is deter- 
mined by solving the Riccati equation and the weighting matrices for the LQR 
problem are similar to those in Ref. 12 except that only two independent 
weighting coefficients are used instead of three. 

In Ref. 21 the weighted sum of the structural mass and control system 
performance index is minimized. 
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Reference 22 treats structural variables as well as coefficients of the 


weighting matrices and orientation of an actuator as design variables. Nu- 
merical results are shown for a two bar truss example. 

In Ref. 23 a nested optimization method is presented for the state feed- 
back control which minimizes the total equivalent mass of the system (struc- 
tural mass plus the mass effect of the control effort). Structural dimensions 
and the coefficient of the control effort are optimized simultaneously to mini- 
mized the objective with a constraint on the mean square of the response. 
Then the control law is determined by solving the Riccati equation with a new 
set of weighting matrices (since the coefficient of the control effort is optimized, 
the performance index is updated for each iteration). 

In Refs. 24 and 25 locations of actuators and sensors are treated in terms 
of (0,1) discrete variables. A utopian multiobjective function containing 
structural mass, control effort and number of actuators is minimized by treat- 
ing structural variables, (0,1) actuator/sensor location variables and open-loop 
gains as independent design variables in optimization. 


1.3 SCOPE OF THE WORK 

In this study the finite element method and linear state feedback are 
combined to formulate the control augmented structural optimization problem. 
A truly simultaneous structural/control optimization scheme is presented in the 
sense that it uses only one set of constraints and one set of design variables 
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(which includes structural sizing, control gain and actuator mass variables). 
As mentioned earlier, this scheme was usually adopted for the output feedback 
control case where the number of gains is relatively small, but by extending 
design variable linking concepts to the control system gains, design space inte- 
gration is achieved for the case of full state feedback control while using a rel- 
atively small number of control system design variables. Several control design 
variable linking schemes are presented and their feasibility and effectiveness 
are shown by solving several examples. 
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Chapter II 

PROBLEM STATEMENTS 

2.1 INTRODUCTION 

In this study a new simultaneous approach to the design of both the 
structural and the control system is presented. The finite element method is 
used to model the structure and linear full state feedback, control is chosen as 
the control law. Dynamic analysis equations from the finite element model 
and the equations of the control system are combined and the design problem 
is formulated as a general nonlinear inequality constrained mathematical pro- 
gramming problem. 

2.2 PROBLEM FORMULATION 

The total mass of the systems has been chosen as the objective function 
and constraints on: (1) dynamic stability (real parts of complex eigenvalues 
or modal damping ratios); (2) damped frequencies (imaginary parts of complex 
eigenvalues); (3) peak transient responses; (4) peak transient control forces; (5) 
control effort; and (6) actutor mass constraints are included in this study. The 
optimization problem treated here can be stated as follows: 

Design a control augmented structural system which minimizes the total 
mass and satisfies various behavior constraints as well as side constraints on 
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the design variables, while treating both structural and control design variables 
simultaneously and independently in the optimization loop. 

In this problem there are three types of design variables. First group 
contains the usual structural design variables (SDV's) such as the cross sec- 
tional dimensions (CSD's). The second group includes nonstructural lumped 
mass design variables such as actuator masses. And Finally the third group 
contains the control design variables (CDV's). Since linear state feedback is 
chosen, the elements of the feedback gain matrix constitute a possible set of 
control design variables. However, this approach is limited because the num- 
ber of control design variables grows very rapidly as the size of the model be- 
comes larger. In order to overcome this limitation several control design 
variable linking schemes are introduced (see Chapter IV). 

The foregoing problem statement can be cast in mathematical form as 
follows: 

Find Y to minimize 

F( Y) (2.1) 

subject to 

Gj{ Y)<0, j= 1 , ... , NCON (2.2) 

with bounds 

\ L { < Y ; < Y-' i = 1 , ... , ND V (2.3) 
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where NDV is the total number of design variables, Y = [Y ( , Y 2 , ... , Y NDV ] T 
is an NDV x 1 design variable vector, F is a scalar objective function, NCON 
is the total number of constraints, (7, is the j-th behavior constraint, and 
Yf, Y y are lower and upper bounds of the i-th design variable, respectively. 
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Chapter III 

' MODELLING 

3.1 FINITE ELEMENT MODEL 

The equations of motion are based on a finite clement formulation. The 
element stiffness and mass matrices for a general frame finite element in local 
coordinates are given in Appendix A. By assembling the element matrices in 
the global coordinates the equations of motion can be written as follows: 

[A/im + [C]{*} + VC]{q] = {F} (3.1) 

where { q } is an ;V x 1 vector of nodal degrees of freedom (DOF), { q } and 
{ q } are first and second time derivatives of { q } respectively, [A/] is an 
.V x N mass matrix, [AT] is an N x N stiffness matrix, [C] is an N x N viscous 
damping matrix, and { F } is an N x 1 load vector. 

It is assumed that the preassigned damping inherent to the structure can 
be represented by a proportional damping matrix which is a linear combina- 
tion of the structural mass and stiffness matrices, i.e., 

[C] = c M [A/] + c k [K] , c A/ , c K constants (3.2) 

There are two kinds of load in the vector {F} of Eq. (3.1): control 
(actuator) forces and external disturbances. With the assumption that the 
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actuator forces and the disturbances act at nodes of the finite element model, 
{ F } can be written as 

{F} = [*]{«} + [>]{/} (3-3) 

where { u } is an M x 1 actuator force vector, M is the number of actuators, 
{/} is an L x 1 vector of external disturbances, L is the number of different 
external disturbances making up a single load condition, and [A] and [e] are 
/V x M and NxL coefficient matrices consisting of the directional cosines 
which respectively relate actuator and disturbance forces to the global coordi- 
nates. 

Now Eq. (3.1) can be written as 

[A/] { q } + [C] { q } + [AT { q } = [A] { « } + [*]{/} (3.4) 

3.2 CONTROL MODEL 

Equation (3.4) can be transformed into the first order state space equation 
{*} = tA 0 Ux} + [B]{u) + [£]{/} (3.5) 

where { x} is a IN x I state vector which is the concatenation of the vector 
of nodal DOF's and its time derivative ({ q }and { q } ), [/1 0 ] is the 2/V x IN 
system open-loop matrix, [5] is the 2N x M system control input matrix, and 
[£] is the 2 N x L system disturbance matrix. The foregoing transformation 
is accomplished by combining the identity { q } = [/]{?} w hh l he result 
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obtained by solving Eq. (3.4) for { q } and then introducing the following no- 
tation: (x} r = , 

r [0] m 1 




-[AfT'c*] -orr'[c] 


(3.6) 


[ 5 ] 


[ 0 ] 

[A/]" 1 [6] 


(3.7) 


[£] 


[ 0 ] 

[A /] _1 l>] 


(3.8) 


where [ ] ' denotes a matrix inverse, and [0] and [/] arc zero and identity 

matrices of appropriate dimensions respectively. 


In control design problems the control law is to be determined. In this 
study linear full state feedback is chosen for the control { u } under the as- 
sumption that all the states (components of { x } ) are available, that is 

{„} = - [//] { x } = - [ [tf p ] [tf v ] ] | R | (3.9) 

where [//] is the M x 2N feedback gain matrix, and [// p ] and [//J are the 
M x N sub matrices containing position and velocity components of [//] re- 
spectively. 


Once { u } and/or [//] are determined, the closed-loop state equation can 
be written as 
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{ -x > = [/(]{x) + [£]{/) 


( 3 . 10 ) 


where the closed-loop system matrix [/l] is 

on - m - m m 

to] in 

- [M] _1 ([/G + [bllHp]) - [M] _1 ([C] 4- [6][// v ]) 


( 3 . 11 ) 
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Chapter IV 

DESIGN VARIABLE LINKING AND INITIAL CONTROLLER DESIGN 
4.1 STRUCTURAL DESIGN VARIABLE LINKING 

In the structural optimization problem some kind of linking scheme is 
commonly used in order to reduce the number of independent design variables. 
In this study two kinds ot cross sections are used for the frame finite elements, 
namely box beam and solid rectangular beam cross sections (see Figure 1). 
For box beam type elements there are 4 cross sectional dimensions which can 
be chosen as design variables, i.e., width (B), depth (H), flange and web 
thicknesses (T2 and T3). For the rectangular beam element there are two 
candidates for design variables: width (B) and depth (H). These cross sectional 
dimensions can be linked in different configurations. In Table 1 possible 
choices for design variables of an element are shown along with the linking 
options within the element. Once design variables for a particular master finite 
element are chosen, it is rather straightforward to link the design variables of 
any other element to those of the master element. 
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4.2 INITIAL STARTUP FEEDBACK GAIN MATRIX 


The purpose of control design variable linking is to keep the number of 
independent control design variables within a tractable range for design opti- 
mization. When any kind of linking scheme is imposed on the feedback gain 
matrix, some design space freedom will be sacrificed and this will usually lead 
to final objective function values that are inferior to those that could theore- 
tically be achieved using a full set of control gain variables with no linking. 
However, the performance of the overall design process will depend not only 
on what kind of linking schemes are used but also on what kind of startup gain 
matrices are used. Three different methods for generating initial startup gain 
matrices, before imposing any kind of control variable linking, are suggested 
here. 

The first initializing method sets the feedback gains arbitrarily (e.g. 
[//] = [0] ) and then carries out a few design iterations without any control 
design variable linking (i.e., all the elements of the feedback gain matrix are 
independent design variables). This allows all the gains as well as the struc- 
tural design variables to change freely for a few iterations in order to find a 
reasonable initial design prior to imposing some linking on the set of M x IN 
control design variables. Even though the unlinked option is used for only a 
few iterations, this can still be a serious restriction, limiting the application of 
the method to small problems. In Chapter VIII this method is only applied to 
small example problems. 
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The second initializing method is to solve the 2N x 2 N nonlinear matrix 
Riccati equation once in order to find the linear optimal control law corre- 
sponding to the initial structural design. The initial gain values obtained from 
the matrix Riccati equation solution are then used to establish fixed ratios be- 
tween the gains that are assumed to hold throughout the design optimization 
process. 

The third approach is the decoupled Riccati equation method, which gives 
an approximate solution to the conventional full order Riccati equation. This 
method uses normal modes to replace the full order Riccati equation by several 
sets of 2 x 2 Riccati equations that have explicit closed form solutions. By 
neglecting the coupling effect in the feedback loop, the gain matrix in the 
modal coordinates can be assembled (with feedback gain vectors corresponding 
to each 2x2 decoupled Riccati equation solution) and transformed to the 
physical coordinates using the normal mode information. 

The second and third initializing methods are described further in the 
following subsections. 

4 . 2.1 Full Order Riccati Equation 

Consider Eqs. (3.4) and (3.5) with the external disturbance terms set to 
zero (i.e. {/} = {0} ): 

C V/] { q } + [C] { q } + IK] { q } = [*] { u } (4.1) 
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{x} = lA 0 -]{x} 4 - [ 5 ]{ w } 


(4.2) 


The optimal control law to minimize a given performance index 

PI = [°°({x} 7 '[£>] { x } + {w} r [fl] {u})dt (4.3) 

J 0 

where [£>] and [/?] are 2N x 2 N positive semi definite and \t x M positive 
definite weighting matrices for states and control forces respectively, can be 
determined from (see Ref. 28) 

{«} = 



= -ixr 1 [Bfcm*} 

where superscript o denotes the initial startup matrix, and the 2 N x 2 N posi- 
tive definite symmetric matrix [Z 5 ] satisfies the following 2N x 2/V nonlinear 
matrix Riccati equation: 

IPJiAA + 1 * J V] + [2] - [/>][«] [«]' l [B] 7 '[/ > ] = [0] (4.5) 

Here { } r and [ ] T represent transposed vectors and matrices, respectively. 

4 . 2.2 Decoupled Riccati Equation Solution 

In this subsection an alternative method which bypasses solution of the 
full order Riccati equation (Eq. (4.5)) is presented. First find the natural fre- 
quencies and normal modes of 
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[A/]{?} + UQ{q} = {0} 


(4.6) 


that is solve the standard eigenproblem 

«, 2 M { Vi } = IK] { v f } i=l ,2, ... ,r (4.7) 

and normalize the modes { v, } so that 

{ v,- } 7 [A/] { Vj } = 5 ip ij = 1 ,2, ... ,r (4.8) 

where r is the number of normal modes retained (r < N) and Sy is the 
Kronecker delta. Let 

{q} = in {z} (4.9) 

where the i-th column of the N x r normal mode matrix [f 7 ] is the i-th normal 
mode { v, } and { z } = [ z u z 2 , ... , z r J r is an r x 1 normal coordinate vector. 
Substituting Eq. (4.9) into Eq. (4.1) and prcmultiplying [f / ] r results in 

{ z } + Diag [c ( ] { z } + Diag [o>, 2 ] { z } = IV] T [6] { u } (4.10) 

where Diag [cj and Diag [<y?] are diagonal rx r matrices whose components 
are c’s and cofs . It is noted that c t = c Sf -f c K coj , in view of the propor- 
tional damping assumption embodied in Eq. (3.2). Equation (4.10) can be 
written in scalar form as follows: 

h + CiZ-i + ajzf = i — 1 ,2, ... ,r (4.11) 

Now assume that the part of the control vector { u } which is related to 
the i-th normal coordinate ( z t and z, ) can be calculated independently and 
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that the resulting { u } is the sum of all the parts corresponding to each scalar 
equation Eq. (4.1 1). 


The foregoing assumption can be stated as follows: 


z i + CjZf + oj] z t = { V,- } T [b ] 


i » > (0 + 



(4.12) 


where 



and 


(4.13) 


(«} = X 5 " 1 * 0 
(=1 

= - [ ( H p ) (l) ( H, } (l) { H p ) (2) { H y } (2) ... { H p } (r) ( H, } W ] 

T 

* [z\ Z\ z 2 h - z r z r 1 

(4.14) 

= - [ { H p } (l) ( H p ) (2 > ... ( H„ } W { } (1) { H „ ) (2) ... { H, } M ] 

T 

X [ z, z 2 ... z { z 2 - ] 

= - [ [Hp [w v ]]r| 

The vector { u }« defined by Eq. (4.13) is an M x 1 control vector which con- 
tains only i-th normal mode information (z h z ( ). Furthermore, { H p 
{ H v }<*> are M x 1 feedback gain vectors which relate { u }« with z h z, respec- 
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tively, while [// p ] and [// v ] denote M x r position and velocity gain matrices 
the columns of which are { H p }<*> and { H v } (,) , respectively. 

In order to recover the initial feedback gain matrix ([//°]) in the original 
coordinates from \_H p ] and [// v ] in the normal coordinates, premultiply Eq. 
(4.9) by [f / ] r [A/] and note that [f 7 ] 7 ” [A/] [ = Q/] » in view of the normal- 
ization imposed by Eq. (4.8), so that 

fz} = m (4.15) 

Then substitute the above equation into the final form of Eq. (4.14) 

{«} = - [// v ](i} 

= - lH p UV] T l MUq) ~ lH v Un T i\ Q{q) (4-16) 

= - [/# ( q } - [/Cl ( R } 

from which it follows that 

[/n = [[//;] [//°] ] (4.i7) 

[/*;] = c/ymTw] (4. is) 

[W°] = [«J['l r [W] (4.19) 

Now the remaining problem consists of finding solutions for the r sets of 
modal gain vectors { H p }d> and { //„}('>. Eq. (4.12) can be transformed into 
the standard first order state space form as 
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{W/} = lAiUw,} + [5,]( { W } (0 + Y {u} 

k^i 


(k) 


0 1 


-°>i ~ c i 


{w/} + 


[ 0 ] 


(4.20) 


( U ) W + £{ U ! 


(k) 


k*i 


where { w t } = |_z, z,J r is the 2 x 1 state vector, [/I,] is the 2 x 2 system 
open-loop matrix and [5,] is the 2 x M system input matrix for the i-th modal 
equation. 


The performance index for the i-th mode ( PI,- ) has the form 

P! i = f°°( { Wj ! T [2/] {>*() + { h } (,Tr [i*,] { uf)dt (4.21) 

J o 

where [£>,] = Diag (0' u , £>22 ) ( Q\ u Q{ 2 > 0 ) is a 2 x 2 weighting matrix for 

the i-th state vector and [7?,] = y, [/] ( [/] is an M x M identity matrix, 
>’; > 0 ) is a weighting matrix for the i-th modal control force vector. Then the 
i-th component of the control { u }W can be determined by 


w 


10 = - [fij-'cs ,-] 7 [/>,]{ IV, ) 


[ to] mV,-! 


_1_ 

y-t 


77 [ [A] 7 ! V; ) p\ 2 

jjin T {y,)Pa*i ~ 77 m T {Vi)Pnii 



(4.22) 
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P\ 1 P\ 2 

where the 2 x 2 positive definite symmetric matrix [R,] = satisfies 

[Pn Pn_ 

the 2 x 2 Riccati equation 

[PRs I/] + [.4,]V/] + [2/3 - iPiiWW~'W T lPii = [0] (4.23) 

Equation (4.23) can be solved in closed form (see Appendix B) and the results 
are 



^11 = c iP\2 + P22 + P\2P22 W i 


where 

W, = (v i ) r M[R i r l M 7 '{»/} = 4-{v,) T [*][ft] 7 '(v i ) (4.25) 

’ l 

By comparing Eq. (4.13) and Eq. (4.22), the i-th feedback gain vectors in 
normal coordinates can be obtained as follows: 

{ H p j (i) = Ibf { V; ( , { ) (0 = [bf { V,. } (4.26) 

By substituting Eq. (4.22) into Eq. (4.12), the i-th closed-loop equation be- 
comes 
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(4.27) 


+ (W + -^■<v j ) 7 '[*]M 7 '{v j }p! 2 ) Z ,. 

= - X ( ^ { Vi } T [6] [ft] 7 { ) (pf 2 z t + p 2 * 2 z*)) 

fc#/ V 

To summarize, the original full order Riccati equation Eq. (4.5) is re- 
placed by r sets of 2 x 2 Riccati equations (Eq. (4.23)) which have explicit 
closed form solutions Eq. (4.24). Then the feedback gain matrix [//] = 
[[// p ] [//„]] in normal coordinates is transformed to [//] in the original co- 
ordinate system by using Eqs. (4.17)-(4.19). 

The method presented in this subsection has a considerable advantage 
over the full order Riccati equation solution approach described in the previ- 
ous subsection. This innovative method is in fact explicit and efficient so that 
this 2x2 Riccati solution procedure can be performed periodically to update 
the fixed ratios initially established by the startup gain matrix. 


4.3 CONTROL DESIGN VARIABLE LINKING (I): ROW-WISE AND 
COLUMN WISE CONTROL DESIGN VARIABLE LINKING 

In this section various linking options for control design variables based 
on row-wise and column-wise linking schemes of the feedback gain matrix are 
presented (see Ref. 26). First the feedback gain matrix [//] can be written in 
various ways as follows: 
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m = [ iH p _ i [//j ] 


L^Ji 

LHJ 2 

L«J,w 


L^Ji 

L«vJi 



1 .Hph 

[H V J 2 



V»p\m 

In Am 



( v 1 

{Up?--- {H p } S 

' {//„}' {«, 


ho a) ■ 


//( 1 ,2 A/) 

H( 2,1) 

H{ 2,2) • 


//(2,2/V) 


//(A/,2) - 


H(M,2N) 


(4.28) 


where [//] is the M x 2/V feedback gain matrix, [// p ] and [// v ] are M x ;V sub 
matrices of [//] containing position and velocity parts respectively, [H\j is the 
j-th row (1 x 2iV) of [//], |_// p J ; and |_//J ; are the j-th rows (1 x ;V) of 

[// p ],[//„] respectively, { H p }' and { H v }' are the i-th columns (M x 1) of 
[// p ] and [// v ] respectively and //(_/,/ ) is the (j,i )-th component of [//] . 


The main ideas underlying the creation of alternative row-wise and 
column-wise control design variable linking schemes are: (1) separation of ve- 
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locity and position parts of the gain matrix; (2) various row and column 
schemes corresponding to actuator and degree of freedom linking; and (3) 
linking schemes based on only allowing changes in various sets of velocity 
gains. Combining the foregoing ideas leads to numerous linking schemes with 
distinct sets and various numbers of independent control system design vari- 
ables (CDV's), ranging from 1 to M x 2 N (see Table 2). 


For example consider option number 5 in Table 2. The feedback gain 
matrix can be written as follows: 


m = 


LtfpJi 

L^vJi 


«iL/£Ji 

«i/+i 1XJ1 

lH p h 

IAJ2 

— 

« 2 lH° p h 

a M+ 2 l H °h 

[ 

■ ' 
s; 




*2 M L<J M 


(4.29) 


Left hand side represents the M x 2N feedback gain matrix in partitioned 
row-wise form ( \_H p \j, [ H v jj represent the j-th rows of [// p ] and [//„], re- 
spectively), and the right hand side has scalar participation coefficients (a/5) 
placed in front of the partitioned rows of the initial startup gain matrix on 
which the linking scheme is imposed (superscript o denotes the initial startup 
matrix). During optimization the a ’s are treated as independent design vari- 
ables (simultaneously with the CSD's) and as they are optimized, the feedback 
gain matrix [//] is optimized in the constrained design space corresponding to 
the fixed ratios established by the rows of the initial startup matrix. 
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4.4 CONTROL DESIGN VARIABLE LINKING (II): 

BLOCK TYPE CONTROL DESIGN VARIABLE LINKING 


In this section a different approach to linking of control design variables 
is introduced and assessed when the decoupled 2x2 Riccati equation solution 
method is used in finding the initial startup gain matrix (see Ref. 27). The 
initial feedback gain matrix obtained by solving r sets of 2 x 2 Riccati 
equations (Eqs. (4.17)-(4.19)) can be rewritten in the following form 


[tf°] = [ [Hp [W°] ] 

= £[ [«;]"’ [h “] (i) ] = 

i= 1 /= 1 


[tfp = £{ H p } (i) { V; } T [.V] = ^[//;] <0 
i= 1 /= 1 


[/C] = £f H, | W { v, } T [W] = 

i= 1 i= 1 


where superscripts (/') indicate that these quantities correspond to the i-th 
Riccati equation. The foregoing equations imply that the [Hffl's and 
[//^P'5 or the [//'’]<*) ' s may be interpreted as basis matrices which can be 
used to generate the initial gain matrix. This suggests that the actual feedback 
gain matrix can be well approximated as a linear combination of these basis 
matrices, namely 


ffl = X! «/ H-33) 


(4.30) 


(4.31) 


(4.32) 
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or 


[>/]=£ [«i [w;] (i) =< i+ ,[w"] w ] < 4 - 34 ) 

/= l 

The whole feedback gain matrix can be linked (Eq. (4.33)), or the position and 
velocity parts of the gain matrix can be block linked separately (Eq. (4.34)). 
During optimization the participation coefficients a 's are treated as independ- 
ent design variables along with the structural sizing variables. It should also 
be noted that these participation coefficients can be further linked with each 
other. 

4.5 PROBLEMS IN CONTROL DESIGN VARIABLE LINKING 

In this section some special aspects of using control design variable linking 
schemes arc discussed. 

When a specific control design variable linking scheme is used, the relative 
values of certain elements in the feedback gain matrix [//] remain frozen 
throughout the design process according to the linking scheme selected. For 
example, if the last option of Table 2 (which links all the elements of [//] to- 
gether) is used, all the ratios among the elements of [//] are invariant. 

But, as mentioned earlier, 2x2 Riccati equations can be re-solved to 
update the fixed ratios between the elements of the gain matrix. This is done 
by finding a new weighting coefficient (y /5 see Eq. (4.21)) for each mode. A 
rational procedure for updating the weighting matrices is described in Appen- 
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dix C. When the new sets of weighting matrices are found, special attention 
is given to preserving continuity of the real part of the closed-loop eigenvalues, 
so that the dynamic behavior remains relatively smooth between the updating 
stages. 


Another problem to be adressed in using control design variable linking 
schemes is that higher modes, which are not directly constrained in the opti- 
mization problem, may become unstable after several design cycles. This effect 
is explained in detail in Appendix D and a way of preventing instability of 
unconstrained higher modes without knowing higher mode information is also 
suggested. 
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Chapter V 

DYNAMIC BEHAVIOR CONSTRAINTS 
5.1 INTRODUCTION 

For the dynamic problems with transient external disturbances the set of 
possible behavior constraints should include; (1) complex eigenvalue con- 
straints; (2) transient response and control force constraints; and (3) control 
effort constraints. 


Every inequality behavior constraint can be written in the form 
g £ S a 


(5.1) 


where g is a measure of a certain behavior and g a is its allowable value. The 
above relation is transformed into a normalized form with respect to its al- 
lowable value such that the constraint function ((?) is always negative. 


where 


G = 4 *> when Za > °’ 

6(2 

G = 1 - -f-, when g a < 0 

6a 


(5.3) 
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5.2 COMPLEX EIGENPROBLEM 


Since the closed-loop [ A~\ matrix (Eq. (3.11)) is not symmetric, 
eigenvalues and eigenvectors are complex and two distinct sets of eigenvectors 
exist for each eigenvalue. For the i-th eigenvalue X, , the right eigenvector 
{ } satisfies 

[/!] {<£,} = X { {</>,•} (5.4) 

and the left hand eigenvector { t /q } satisfies 

!^ i ) r M = (5.5) 


and 

Xj = °i + w di (5.6) 


where cr, and a> di are real and imaginary parts of the i-th complex eigenvalue. 
Also the modal damping factor ^ of the i-th mode is defined as 


^ = 



+ 



(5.7) 


These two sets of right and left hand eigenvectors are normalized such that (see 
Ref. 29) 

(<M r {0/} = 4 (5.8) 


and 

{ tj ) T S,J (5.9) 
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where d t is a normalizing scalar constant for the i-th right eigenvector and b tj 
is a Kronecker delta. After solving both right and left eigenproblem, and 
normalizing the eigenvectors, it can be shown that the following matrix re- 
lations represent valid identities: 


[.d] [d>] = [O] [A] 

m 7 [d] =[A][^] 7 

m 7 [o] = [/] 

[V] 7 [d][cD] = [A] 


(5.10) 


where [<D] = [ <£, (f) 2 ... 4> 2R ~\ is a 2Nx2R eigenmatrix in which the i-th 
column is the i-th right hand eigenvector, [4 / ] = - 'I'ir ] is a 

2N x 2 R eigenmatrix in which the i-th column is the i-th left hand eigenvector, 
[A] is a 2 Rx 2 R diagonal matrix which contains eigenvalues on its diagonal 
and 2 R is the number of eigenvalues and eigenvectors considered ( R < AO- 


5.3 DYNAMIC TRANSIENT RESPONSE ANALYSIS 

Time dependent response and control force constraints are replaced with 
a finite number of peak constraints by first finding the corresponding peak 
times using the adaptive one dimensional search method described in Ref. 30. 

In order to calculate the transient response of the given system Eq. (3.10) 
should be integrated in the time domain. Since Eq. (3.10) is coupled and usu- 
ally large, a set of complex eigenvectors is chosen as a basis to diagonalize and 
reduce the dimensions of the original equation. Let 
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{•* } = [®] { 1 i 


(5.11) 


where [<t>] is a 2/V x 2R ( R < N ) right hand eigenmatrix and { rj } is a 
2Rx 1 complex normal coordinate vector. Substituting Eq. (5.11) into Eq. 
(3.10) and premultiplying by ( 2Rx2N left hand eigenmatrix trans- 

posed) yields 

m 7 [<t>] {v } = 4- [V] 7 [£]{/} (5.12) 

and introducing the identities in Eq. (5.10) leads to 

il! - [A] { I) } + m 7 ' [£](/[ (5.13) 

Equation (5.13) is equivalent to the following sets of 2 R first order scalar dif- 
ferential equations 

Vi = 2 i rj i + {( A/} 7 [£]{/}, /= 1,2, ... ,2 R (5.14) 

Integrating Eq. (5.14) with respect to time gives 


V iiO = 


to) V iVo) 


+ 




A,-} T [£] {/(T)}dr 


(5.15) 


where t 0 is an initial time and t is a specific time of interest. Here {/( r) } is 
assumed to be expressed in terms of a truncated Fourier series and 
polynomials over a specified period of time i.e., 
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(5.16) 


{At)} = j]({ FC] k cos n k t+ {FS) k s\ n n t t)+Y J { Fp }p ,P ’ 

k = 1 p=0 

when 0 < t < tf 

{/( 0} = o . when * > 7 

where tj- is the specified time interval during which disturbances arc applied, 
N q is the number of different driving frequencies, {FC} k and {FS} k are 
L x 1 vectors which contain cosine and sine components corresponding to the 
k-th driving frequency Q k and N P denotes the highest order polynomial term 
considered and the Lx 1 vector { FP} p corresponds to the p-th order 
polynomial. Then substituting Eq. (5.16) into Eq. (5.15) yields 


0) 


/ 


+ 


{« A/} r [E] f e k ‘ {t T) V ({7C} fc cos a k x+{FS} k s\n Cl k x) \ dx 

J ° v= i ’ 


N P 


(5.17) 


+ { } T IF] ^ {FP} p X p ) dx when 0 < t < t f 


P= 0 


rn(t) = 


when t > t 


7 


More detailed expressions for the rj's are given in Appendix E. After calcu- 
lating the rj-s, i= 1, ... ,2 R the state vector (.v(f)} can be recovered from 
Eq. (5.1 1) 

{ jt } = [<t>]{')} U") 


and the control force vector { u (t) } (see Eq. (3.9)) is given by 


34 


{ u } = -[//]{ x } = - [//] m { v } 


(5.18) 


Also the transient acceleration can be calculated by differentiating Eq. (5.11) 
with respect to time and using the relation Eq. (5.13) to give: 

{*) = mUl) = [4>] ([A] { T] } + [H'] 7 ’ [£]{/}) (5.19) 


5.4 CONTROL EFFORT 


The control effort can be defined as 


Control Effort 



roo _ (5.20) 

{u } 1 [R CE ]{u}dt 


where [ R CE ] is an M x M positive definite symmetric weighting matrix. 


When the time tj (during which the external disturbance force is applied) 
is small, the first term of the right hand side is negligible compared with the 
second term. So in this work the control effort (CE) is defined as 


CE = f°°{ U } T [/? CE ]! u }rf, 

J '/ 

= p{.v}' r [ff] r [tf CF ] [//]{*}<* 

J '/ 

= r { x ) T IQceI { * } * 


where [<2 CE ] = C^] r [^cE] OT- 


(5.21) 
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Consider the closed-loop equation (Eq. (3.10)) when t > t f 


{ x } = [/l] {x } , t > t f 


(5.22) 


Substituting Eq. (5.11) into Eq. (5.22), premultiplying by [ v F] r , and using the 
identities in Eq. (5.10) leads to the following expression 

{</} = [A] {if}. ' £ '/ < 5 - 23 ) 

Substituting Eq. (5.1 1) into Eq. (5.21) above gives 

CE = f°°{ n } r [<i>] r [a CE ] [4>]{>l} dt 

'/ (5.24) 

poo ^ ~ 

{ri} T [Q-]{ri}dt 

Jl f 

where [0] = [ c b] r [(2ct] [$] ■ From Eq. (5.23), the complex response 
{ rj (0 } can be written as 


(!!(<)) = (!)((/)). 


' * '/ 


Substituting Eq. (5.25) into Eq. (5.24) yields 


ce = {,( V )) r rvtf* 

J tr 


[A](/-r 


= r e CA fU-V&eW-Vdt 

•* w 
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[A ] 7 [IT] + [H>] [A] 

= f 00 [A] 7 '^ r ^^/) [ g ] ,[A](r-/ / )^ 

’/ 

+ f°V A l r ( f “ l f ] [0] e [A] {t ~ t f ) [A] dt 

J 'f 

= * (5.28) 

CO 

'/ 

= [0] - W 

= - [5] 

>r 

[A] [W~\ + [#][ A] = - [£] (5.29) 


which is a Lyapunov equation. Therefore the control effort of Eq. (5.26) is 
given by 


CE = {n(lf)) T V ihiriOf)} (5.30) 

where [1L] satisfies Eq. (5.29). Equation (5.29) can be solved clement by ele- 
ment since [A] is a diagonal matrix and \_Q ] is a symmetric matrix, i.e., 



Qu 


Xi + Xj 


(5.31) 
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where W tj and Q i} are the ( i,j)-th element of [H / ] and [g] ( Q tj - 
{ ■(j) i } r [0 CE ] { 4>J }) and ( 4>t } is the i' th complex right hand eigenvector. 


5.5 ACTUATOR MASS CONSTRAINTS 

The actuators, located at specified structural nodes, can be sized so that 
they produce control forces or torques required to suppress the vibration. In 
Ref. 31, it is assumed that actuator masses are fixed during one cycle of opti- 
mization, and they are updated after finding peak control force or torque val- 
ues of the new design according to functional or empirical relations between 
the maximum peak control forces or torques and the required actuator mass. 

When the closed form solution for the transient peak control forces or 
torques is available, the relation between the actuator mass and the control 
forces or torques can be mathematically stated in a constraint form as follows: 

m., a c, | lit) |% 0<(</ m „ (5.32) 

or 

m A > c, | u{t-)\ c \ j= I , ... , NPEAK (5.33) 

where m A is the mass of the actuator, / max is the time interval of interest, tj is 
the j-tli peak time for control forces or torques, NPEAK is the number of 
control force peak times, and c, , c 2 are constants relating the peak control 
forces and required actuator masses. Eq. (5.33) can be normalized such that 
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1 4tJ) f 2 

~A 


- 1 < 0 , 


j= l,... ,N PEAK 


( 5 . 34 ) 


Chapter VI 

SENSITIVITY ANALYSIS 
6.1 INTRODUCTION 

In order to generate approximate problems first order system sensitivity 
information with respect to both structural and control design variables is re- 
quired. Since it is assumed that the external loads are expressed in terms of a 
truncated Fourier series and polynomials over a specified period of time, it is 
possible to calculate all of the (first order) behavior sensitivity derivatives an- 
alytically. 


6.2 FUNDAMENTAL SYSTEM MATRIX DERIVATIVES 


Derivatives of [//] with respect to control design variables are obtained 
directly and derivatives of [A/1, [C] and [ K] with respect to structural vari- 
ables are obtained analytically from the finite element formulation. From this 
information analytic sensitivities of the system [/<] and [£] matrices (see Eq. 
(3.10)) can be obtained with respect to an arbitrary intermediate or direct de- 
sign variable a , as follows: 


[ 0 ] [ 0 ] 

ML = 

da o[A 2 i~i d[A 22 ~\ 

da da 


( 6 . 1 ) 
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and 


£[g] 

5a 


[0] 

_ da J 


where in general 


(6.2) 


5 [^ 2 |] 

da 


ca 


d [A/] 


-1 


5a 


-(ra + [*][«„]) - on 


<5[/y 




(6.3) 




5[£ 2 ] 

5a 


c[M] _l 


c<a 


w 


(6.4) 


It should be noted that the matrices [5] and [e] are constant with respect to 
any design variable and that according to the type of the design variable a , 
structural or control, some matrices are independent of that kind of design 
variable. In other words, the derivatives of [A/], [C] and [/fj with respect to 
the control design variables are zero and those of [// p ] and [//„] with respect 
to the structural design variables are also zero. The derivative d\_.\f\~ x lda > s 
calculated by differentiating [A/] [A/] -1 = [/] = constant , which results in 


aoo~' 

da 


- [A/] - 


5[A/] 

da 


[A/]- 1 


(6.5) 
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6.3 EIGENPROBLEM SENSITIVITIES 


Eigenvalue sensitivities are needed not only because eigenvalues them- 
selves can be behavior constraints but also because they are required to calcu- 
late other response sensitivities (see Eq. (6.18)). Furthermore, in order to 
obtain precise transient response sensitivities, eigenvector sensitivities are also 
calculated using an analytic method (Refs. 32 and 33). 


Differentiate Eq. (5.4) and (5.5) with respect to a and rearrange terms 


(M - X, [/]) 


da 


da 


dX; 




(6.6 .a) 




T 


da 


([-4] - *,[/]) = - (<M 


t ML - B.f*,} 7 ' 

da da 


(6.6 .b) 


or 


[ZJ 


da 


- {?>} 


) 


T 


o a 


CZJ = - {G t } 


T 


where 


[zj = [ a ] - ;.,[/] 


or-ll <M; 


( ^ -,T (I ,7 5[/I] f , ,7 

< c -> = 




(6.7.6) 


( 6 . 8 ) 


(6-9. a) 


d a 


(6.9.6) 


Premultiplying Eq. (6.6. a) by { \]/ t } r and using Eq. (5.5) and Eq. (5.9) leads to 
the following expression for eigenvalue derivatives 


dX: do-. do)j: T d[A~\ 


da da da 


or 


do : 


oa 


f >7 ADO , x y ( , y 7 
\VifR — t Vi fR ~ ( }/ ~ ' a ■ { <f>i }[ 


da 


and 


do} di f , ,T d[/4] 


7 b\.A] 


doc 


= i+ i VR- J £ L i4>ih+ ihVi^itih 


doc 


(6.10) 


( 6 . 11 ) 


( 6 . 12 ) 


where subscripts R and / represent real and imaginary parts. The sensitivity 
of the modal damping factor is calculated by differentiating Eq. (5.7) with 
respect to a as follows: 




m d\ °i 


da} di 

da 


do i 
da 



da 


f 2 

.of + 



3/2 


(6.13) 


The derivative of the i-th eigenvector can be expressed as the sum of ho- 
mogeneous and particular parts, as follows: 


d_W 

doc 


«,{*/} + {Vi) 


(6.14.a) 


H+i} T 

doc 


= bii+tf + {Uf 


(6.14.6) 
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The particular solutions { V t } and { U t } T satisfy 
[Zj]{ K,( = - {Fj} (6:1 5. a ) 

{ V;} T {Z^ = - (Ci) 7 (6.15.6) 

where [ZJ is matrix [Z,] with the r-th column and row replaced by { e , } and 
{e r } T respectively, { e r } is a vector containing 1 in the r-th element and 0's 
elsewhere, { F t } and { G , } are vectors { F t } and { G t } respectively with the 
r-th element replaced by zero, where r is the location of the maximum absolute 
element in { <f >, } . The coefficients of the homogeneous parts, a { and 6„ can 
be calculated by differentiating the normalization conditions Eq. (5.8) and Eq. 
(5.9), noting that d ( is a constant and i = j as follows: 

j^{h) T W) = o. ±{ {tl ) T {*,}) - o 

hence (in view of Eq. (6.14)) 

{<h) T (<•,{<!>,) + W) - o. 

(&,■{ <h) T + { {<#»,} + l*,) T {*,",} + ! Ki) = 0 


so that 


( 0/ ) T l<h) d ' 


( 6 . 16 .fl) 


6 ,= - a, - {<h) T { v i) 


( 6 . 16 . 6 ) 
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6.4 TRANSIENT RESPONSE SENSITIVITIES 


For an arbitrary peak time t, peak response sensitivities can be calculated 
by differentiating Eq. (5.1 1) as 


d {-*(0} 

da 




a[»] 

da 


+ [O] 

da 


(6.17) 


In the above equation, the only term not previously derived is d{rj (t) }/da and 
it can be obtained by differentiating Eq. (5.15) with respect to the design var- 
iable a, assuming that the external disturbances are independent of the design 
variables, i.e. Q k , {FC} k , ( FS } k , k = 1, N a and {FP} py p = 1, N P are constant 
as follows: 


£M0 

da 


JL J <o) 

8a 


v iU„) + JV‘ ( ' ,) !^} r [£]{/(T)}<* 


+ J^('- V' ( '~ ,) {'M 7 '[£] {/(*)} <#r 

+ r^(r-T) 

( o 


(6.18) 


da 


[£] {/(T)}dT 




Once d{rj (t)}/da is calculated, sensitivities of control force and acceleration 
can be obtained by differenting Eqs. (5.18) and (5.19), respectively. 
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A more detailed derivation of <fy ( (r)/da is given in Appendix E. It should 
be noted that the above calculation is carried out only at previously identified 
peak response or peak control force times. 


6.5 CONTROL EFFORT DERIVATIVE 


The sensitivity of the control effort with respect to the design variable a 
is obtained by differentiating Eq. (5.30), 


d(CE) 

da 




dot 


In the above expression d{rj( tf)}\d a can be calculated according to the pro- 
cedure presented in the previous section for / = t f and the only unknown term 
is d[lF]/da and it can be determined by differentiating Eq. (5.31) as follows: 


a fry 

da 


- M ' ■) 

da \ Aj + Aj J 



1 ^Qjj 
Aj + Aj da 


(h + ijr 


dAj 

da 


+ 


dAj \ ~ j dQjj 

) Q,J ~ (2, + Aj) da 


but 


a, = {0,} r cecE]{^> 


( 6 . 20 ) 


( 6 . 21 ) 


so that 
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( 6 . 22 ) 


d Qj ( x , T d tQcE~] f ± , 

= { </>/} T {<£/} 


5a 


+ 


5a 

T T ^{0/} 
^ — CQcrM^/> + {<M 7 [2 CE ]-^ L - 


3 { 0 ,} 


6.6 ACTUATOR MASS CONSTRAINT DERIVATIVE 

The sensitivity of the actuator mass constraint with respect to an arbitrary 
design variable a is obtained by differentiating Eq. (5.34) so that 


dG 5 
= C 1 


«(»; ) l' 2 


da 


da 


m x 


(6.23) 


which gives 


dG 

da 


-- c. 


m. 


~ c 2 




c 2 -l 


5«(r.- ) 


5a 


+ 


K'l" £Ur) 


(6.24) 


when w(r‘) > 0 and 


dG 

da 


= c. 




~^7 Cl 


(6.25) 


when u(tj)<0. In the above expression, 


Mt}) 

5a 


is the sensitivity of the peak 


control forces, and * s n °t>zero only when a is the corresponding 

actuator mass variable. 
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Chapter VII 
OPTIMIZATION 


7.1 INTRODUCTION 

Various approximation concepts such as structural and control design 
variable linking presented in Chapter IV, temporary constraint deletion and 
intermediate design variables (Refs. 34 and 35) are used to replace the original 
optimization problem stated in Chapter II by a series of explicit approximate 
problems. With the first order sensitivity information derived in Chapter VI, 
linear, reciprocal or hybrid approximations can be made with respect to either 
direct or intermediate design variables, even though the approximate design 
optimization problems are always solved in an integrated design space that 
spans the actual structural CSD's and the participation coefficients of the 
linked control gains. Each approximate optimization problem has its own 
lower and upper bounds on the design variables determined by given move 
limits and retains a set of constraints which are active and potentially active 
for the approximate problem. 
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7.2 INTERMEDIATE DESIGN VARIABLES 


In order to generate a high quality approximate problem intermediate 
design variables are used. For frame elements it is known that the section 
properties A, I y , I 2 and J are a good choice for intermediate design variables. 
This follows from the fact that the elements of the stiffness and the mass ma- 
trices are linear functions of these section properties (see Appendix A). 

Control design variables are used directly in the generation of the ap- 
proximate problem because the system matrices are linear functions of the 
gains. After some numerical experimentation, most behavior constraints used 
in this work are found to be adequately approximated by the linear approxi- 
mation with respect to control design variables. So when generating the ap- 
proximate functions the linear approximation is always chosen for control 
design variables although the hybrid approximation can be used with respect 
to intermediate structural design variables. 


7.3 TEMPORARY CONSTRAINT DELETION 

In order to Find a reduced set of constraints for each approximate problem 
the following rules are used: (1) all constraints which are greater than a given 
cutoff parameter are retained; and (2) when all the constraints in one specific 
category, for example when all the peak actuator force constraints are less than 
the given cutoff parameter, the one which has maximum value among them is 
retained so that the effect of this specific category will not be neglected. 
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After determining which constraints are retained for the approximate 
problem, sensitivities of those retained constraints along with the objective 
function are calculated analytically as shown in Chapter VI. 


7.4 APPROXIMATE PROBLEM 

With the information acquired from the analysis and sensitivity analysis 
phase each approximate optimization problem can be formulated as follows: 

Find a design variable vector Y to minimize 

F ( X ( Y ) ) (7.1) 

subject to 

G,(X(Y))<0, J e Q r (7.2) 

with bounds 

Yf < Yj < Y f i= 1 ,...,NDV (7.3) 

where NDV is the total number of design variables, Y is a vector of design 
variables {NDV x 1) , NIDV is the total number of intermediate design vari- 
ables, X ( Y ) is a vector of intermediate design variables ( NIDV x 1) , F (•) is 
an approximate objective function, G } (•) is the j-th approximate constraint, 
Yf- and Yf 7 are lower and upper bounds of the i-th design variable which are 
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determined by a given move limit, and Q R is the retained set of constraints for 
the approximate problem. 

Either linear , reciprocal or hybrid approximation can be used for F (•) 
and Gj (•). The details of the hybrid approximation follow. Let / (•) be any 
approximate function ( F (•) or Gj (•) ), then during the approximate opti- 
mization cycle / (•) is 


/ (X ( Y )) = f(X ( Y 0 )) + 


NIDV 

£• 




3/(X ( Y 0 )) 

ax* 


B k (X( V)) 


- Jo + 


NIDV 


£ C ot S A (X(Y)) 
k=l 


(7.4) 


where Y 0 , f 0 =/(X ( Y 0 )), C ok = df{X( Y 0 ))fdX k , k= I , ... , NID V are values 
at the beginning of the approximate problem which remain constant during the 
cycle, and 


B k {X(Y)) = X* ( Y ) - X*(Y 0 ), 
**<X( Y))= - (X* ( Y„ )) 2 U 1 


when C ok > 0 


(7.5) 


X* ( Y ) X*(Y 0 ) 


, when C ok < 0 


and derivative of / (•) with respect to the i-th design variable Y ( - is 


~ XlDV 

cf (X ( Y )) v V dB k(X(Y)) 

■ = /j C ok - 


C Y: 


k= 1 


dY: 


(7.6) 


where 
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when C ok > 0 


(7.7) 


as*(X( Y)) 3X*(Y) 

«v , OY, ’ 

3B*(X(Y)) /X t (Y„) 

3Y, V X *(Y> 

Note that when no 
N ID V = ND V and X ( Y ) = 


2 «t(Y) 



when 


intermediate design 
Y or X; = Y f , i= 1,2, 


Cok < 0 

variables 
... , ND V. 


are 


used, 
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Chapter VIII 

NUMERICAL RESULTS 


8.1 INTRODUCTION 

The control augmented structural optimization solution method described 
in the previous chapters has been implemented on the IBM 3090 main frame 
computer at UCLA. Numerical results for several examples, which illustrate 
the effectiveness of alternative control design variable linking schemes, are 
presented here. In most example problems, CONMIN (Ref. 36) is used as the 
optimizer, and it is conservatively assumed that the passive damping is zero 
( c M = c K — 0, see Eq. (3.2)), unless otherwise specified. The convergence cri- 
terion used in these examples is that the relative change in the objective func- 
tion values between two sets of consecutive design iterations should be less 
than 0.1 percent. 


8.2 EXAMPLE I - CANTILEVER BEAM 

The first example is a cantilever beam as shown in Figure 2 free to deflect 
in plane (E =7.1 x 10 6 N /cm 1 , p = 2.768 x 10 3 kgjcm 3 ). It has a box beam 
type cross section and is modelled by 10 equal length (100 cm) finite elements 
resulting in 20 degrees of freedom (N = 20). A concentrated mass of 200 kg 
is located at the middle of the beam and a single translational actuator 


53 


( M = 1 ) having a mass of 4 kg applies control force u(t) at the tip. A tran- 
sient half sine pulse (with a magnitude of 4000 N at a frequency of 3.9 Hz) is 
applied at the tip translational degree of freedom for a time interval 0 < t < 
0.1282 seconds. The dimension of the first order state equation (see Eq. (3.10)) 
is 40 (2 AO, so the total number of complex modes is 40. Transient responses 
arc considered for the time interval 0 < t < 2 seconds and the lowest 20 out 
of the 40 complex modes (according to the absolute value of the imaginary 
parts of the complex eigenvalues) are used to calculate the transient responses. 
The normalizing constants for the right hand complex eigenvectors (d/s, see 
Eq. (5.8)) are chosen to be 10 5 . 

For each finite element the width (B) and depth (H) are fixed (B = H = 40 
an), but the flange and web thicknesses are free to change (subject to side 
constraints 0.5 cm < 12, T3 < 10.0 cm) resulting in two design variables per 
finite clement. Initial thicknesses are all set to 5.0 cm. 

Total mass is taken as the objective function and it includes: (1) the fixed 
mass at midspan (200 kg); (2) the fixed actuator mass at the tip (4 kg); and (3) 
the variable structural mass. All the finite elements are linked resulting in 2 
structural design variables (SDV's). Hybrid approximations in terms of the 
structural intermediate design variables (sectional properties) are used while 
the approximations with respect to the control design variables (CDV s) are 
linear. Behavior constraints are imposed on: (1) the real part of all the re- 
tained complex modes ( a t < —0.5 , i = 1,... 10): (2) the lowest damped fre- 
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quency (oj di > 4.0 Hz ); (3) the peak tip displacement, 

(|< 7 19 (t)| ^ 10.0cm); and (4) the peak actuator force ( | u(t) | < 1000 N). 

8.2.1 Control Design Variable Linking (I) 

In this section, row-wise and column-wise control design variable linking 
schemes which are described in section 4.3 arc applied to the cantilever struc- 
ture. The maximum possible number of control design variables is relatively 
small (.1/ x 2 N = 40), so all available control design variable linking options 
are tried (see Table 2). Since there is only one actuator, there are 5 distinct 
control design variable linking options for the same problem instead of 10 (for 
M = 1 each column of the gain matrix contains only one element, so the 
options 1, 2, 5, 7 and 6 are identical to 3, 4, 8, 9 and 10, respectively, see Table 
2 ). 


Initial startup feedback gains are computed by solving 10 sets of 2 x 2 
Riccati equations corresponding to the lowest 10 normal modes (r = 10). The 
control weighting coefficients y-s are set to 1/300 and the 2x2 state weighting 
matrices are chosen to be [£>,] = Diag ( o}, 1 ) / = 1 , ... , r so that the first term 
of Eq. (4.21) represents a total (strain and kinetic) modal energy. 

Five distinct control design variable linking options are imposed on the 
startup gain matrix from the beginning. Move limits of 50 to 60 percent for 
structural variables and 90 percent for control variables are used. Iteration 


histories and final structural designs are given in Tables 3 and 4 respectively, 
and are also displayed in Figures 3 and 4. 

In case I the feedback gain matrix is totally unlinked (or equivalently 
columns of [77] are linked, Option 1 or 3 of Table 2). The number of control 
design variables (CDV's) is 40. This case allows complete freedom in the 
control design variable space and the final design is a near minimum gauge 
design with respect to the cross sectional dimensions (see Table 4). In case 2 
the position part of the feedback gain matrix is fixed and the elements of ve- 
locity part are chosen as design variables (or columns of velocity part are 
linked, Option 2 or 4 of Table 2), resulting in 20 control design variables. In 
this case the convergence is relatively slow and the optimum mass is high con- 
sidering the freedom given to the design space. This is because even though 
the elements of the velocity part of the feedback gain matrix can change freely, 
the fixed initial position gains are not appropriate. This can be attributed to 
the fact that when the position gains are fixed at their initial values, the peak 
displacement and frequency constraints must be satisfied by increasing the 
structural stiffness (once the peak actuator force constraint has become active). 
In case 3 the position and velocity parts (or the rows of the position and ve- 
locity parts) of the feedback gain matrix are linked (Option 5 or 8 of Table 2), 
which leaves 2 independent control design variables. The results for this op- 
tion are remarkable because it takes only 10 analyses to completely converge 
to a near minimum gauge design using only 4 independent design variables (2 
SDV's and 2 CDV's). In case 4 the position part of the feedback gain matrix 
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is fixed and the velocity part (or the rows of the velocity part) of the feedback 
gain matrix is linked (Option 7 or 9 of Table 2) which leaves a single control 
design variable. In case 5 the row of the feedback gain matrix (or the the en- 
tire feedback gain matrix) is linked (Option 6 or 10 of Table 2) which leaves 
also only one control design variable. Cases 4 and 5 exhibit similar conver- 
gence histories and final designs except for the position part of the feedback 
gain matrix (which is frozen in case 4). As one would expect when seeking a 
minimum mass design, the final web thicknesses in all cases are at the lower 
bound value (i.e. T3 = 0.5 cm , see Table 4 and Figure 4). 

For comparison purposes the results for cases 1-5 are summarized in Ta- 
ble 5 along with the active constraints at the final design. Most cases con- 
verged smoothly and the maximum difference in the final objective function 
values achieved was as less than 15 percent. As expected the final design mass 
becomes larger as more restrictive linking is imposed, but the convergence be- 
comes more robust and the total number of analyses required to obtain con- 
vergence is reduced (e.g., see cases 3-5). 

8.2.2 Different Initial Gains for Generating Startup Gain Matrix 

In this section different arbitrary initial feedback gain matrices are used 
to find startup gain matrices in order to examine the first initializing method 
(see section 4.2). It will be recalled that this method chooses arbitrary initial 
control gains to start and then allows a few totally unlinked iterations before 
imposing any control design variable linking. Four different initial feedback 
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gain arrangements are chosen and 3 iterations are allowed without any control 
design variable linking (40 independent CDV's, option 1 or 3 of Table 2). 
Then the position and velocity parts of the gain matrix are linked (2 inde- 
pendent CDV's, option 5 or 8 of Table 2). 

In case 6 the initial gain matrix is null ([//] = [0]). In case 7 the initial 
gain matrix is chosen in a manner similar to that which would be used if direct 
output feedback control was being employed (//(1,19) = 5 kg! sec. 2 , 

H{ 1, 39)= 100 kgjsec., H{ 1, /)= 0, elsewhere). In cases 6 and 7, in addition to 
the larger move limits for the control design variables (i.e., 190 percent), abso- 
lute move limits are used for small elements (namely, for gain elements less 
than 1 move limits are set to 100) during the unlinked iterations. In case 8 the 
initial gain matrix is calculated by solving the full order Riccati equation with 
weighting matrices [(>] = Diag [ [£], [A/] and [/?] = ^ [/] . In case 9 

the initial feedback gain matrix is obtained by solving 10 sets of 2 x 2 Riccati 
equations followed by three unlinked iterations to generate the startup matrix. 

Iteration histories and final structural dimensions are given in Tables 6-7 
and Figures 5-6. All four cases converge to similar final mass values. For 
comparison the results for cases 6-9 are summarized in Table 8 along with the 
active constraints at the final design. In this example the final web thicknesses 
in all cases take on their lower bound value (i.e. T3 = 0.5 cm, see Figure 6). 

The results for cases 6-9 show that the initial gain matrix used to generate 
a startup gain matrix (via three iterations without any CDV linking) does not 
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have a significant effect on the Final structural design or the minimum mass 
values achieved. It is important to note that the method used in this section 
can not be extended to larger problems because as the number of CDV's in- 
volved increases (i.e., according to number of CDV's = M x 2 N ), the three 
totally unlinked optimization cycles required to generate the startup gain ma- 
trix become intractable. 


8.3 EXAMPLE 2 - ACOSS FOUR STRUCTURE: CDV LINKING (I) 

The second example is the ACOSS FOUR structure shown in Figure 7. 
Several studies have been made on this model (see Refs. 11, 13-16, 18 and 19) 
since it represents one of the simplest configurations for a 3 dimensional space 
structure. It consists of twelve truss elements and has 12 degrees of freedom. 
Nondimensionalized Young's modulus and mass density are 1 and 0.001 units, 
respectively. The edges of this tetrahedral truss are 10 units long, and the six 
supporting links are 2.83 units long. An actuator is located in each of the six 
supporting links (elements 7 through 12) introducing control forces that act 
along the supporting links, and 4 nonstructural masses of 2 units each are at- 
tached at nodes 1-4. The nodal coordinates are given in Table 9 (from Ref. 
11 ). 


The design objective is to minimize the structural mass satisfying given 
behavior constraints on the closed-loop eigenvalues. Truss member areas are 
the design variables and no structural linking scheme is used, therefore, the 
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number of independent structural design variables is 12. Initial truss member 
areas are given in Table 10 along with the final areas for all 4 cases. Lower 
and upper limits on the areas are 10 and 1500, respectively. The most restric- 
tive control design variable linking scheme (option number 10 of Table 2) is 
chosen which leaves only one independent control design variable (elements of 
the feedback gain matrix are all linked). Hybrid approximation in terms of 
areas and linear approximation in terms of the control design variable are used 
to construct a sequence of explicit approximate problems. DOT (Ref. 37) is 
used as the optimizer. 

Initially two control feedback gain matrices (Initial Design A and B) are 
found by solving lull order and 2x2 Riccati equations, respectively. Initial 
design A is found by solving the full order Riccati equations with identity 
weighting matrices for both states and control forces (El23 C^l 2 /V> 

[/?] = [/] v „ see Eq. (4.3)). In Tables 1 1 and 12 initial feedback gains, natural 
frequencies and closed-loop eigenvalues of the initial design A are given. Ini- 
tial design B is obtained by solving 12 sets of 2 x 2 decoupled Riccati equations 
for the initial control gains with identity control weighting matrix ( 
[/?.] = [/] w or y’s = 1), and diagonal 2x2 state weighting matrices ([&] = 
Diag ( o of, 1), < = 1, ... , 12 , see Eq. (4.21)). The normalizing constants d's for 
the right hand eigenvectors (see Eq. (5.8)) are chosen to be 1. Again in Tables 
13 and 14 initial feedback gains, natural frequencies and closed-loop 
eigenvalues are shown. In the following subsections 4 cases are investigated 
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where different sets of behavior constraints and initial startup control gains are 
used. 

8.3.1 Inequality Constraints 

In the first two cases, inequality behavior constraints are imposed on the 
first damping ratio (<!;, > 0.15), and first two damped frequencies (gj^, > 1.341 
and ( 0^2 > 1.6, see Ref. 19). Cases 1 and 2 start from the same member areas 
but different control gains, namely, from initial design A and B, respectively. 
Move limits used for the structural variables are fixed at 50 percent, and for 
the control variable the move limit is reduced by 70 percent after each iteration 
starting from 70 percent (case 1) and 50 percent (case 2). 

Final member areas as well as the final value of the single control variable 
arc given in Table 10. The value of the control design variable represents the 
ratio of the final feedback gain elements with respect to the initial gain ele- 
ments, so the actual final feedback gain matrix is the initial feedback gain 
matrix multiplied by the final control design variable in each case. Natural 
frequencies and closed-loop eigenvalues of the final designs are given in Tables 
15 and 16. In both cases the two damped frequency constraints are active at 
the final design but the first modal damping factor is not critical (i.e., case 1, 
0.16465 >0.15 and case 2, 0.16745 > 0.15). 

Iteration histories of cases I and 2 are similar to each other (see Table 19 
and Figure 8). Total design masses include 8 units for the fixed nonstructural 
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masses as well as the structural member masses. Final truss areas are shown 
in Figure 10. 

Cases 1 and 2 show similar area distributions and final mass values (see 
Figure 10) even though the final control gains are different. The final struc- 
tural masses obtained here (14.52 and 14.12, see Table 10) are lower than the 
result reported in Ref. 19 (18.56) by more than 20 percent. This reduction in 
the final mass can be attributed primarily to integration of the structural and 
control design variables, even though the most restrictive control variable 
linking option is used. It may be noted that less restrictive control design 
variable linking options do not reduce the final mass much further in this ex- 
ample. 

8.3.2 Equality Constraints 

The next two cases (cases 3 and 4) are similar to cases 1 and 2 but have 
different behavior constraints. One inequality constraint and five equality 
constraints are imposed^ the first 4 modes arc constrained to ha\e the same 
damping ratios £,= 0.1093, / = 1 , 4 ; and the first two damped frequencies 
are constrained as follows, co d] = 1.34 and > 1.5 (see Refs. 18 and 15). 

Each equality constraint is replaced by 2 inequality constraints that define 
a small interval, namely 0.1093 < £, < 0.11,/= I , — ,4 and 
1.34 < co d] < 1-345 , resulting in 11 inequality constraints. Cases 3 and 4 
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start from initial design A and B, respectively, and both structural and control 
design variables have 70 percent move limits. 

Final designs, natural frequencies, closed-loop eigenvalues are given in 
Tables 10, 17 and 18. All the constraints are satisfied with less than 0.5 per- 
cent constraint violation. 

Iteration histories are shown in Table 19 and also in Figure 9. Final truss 
areas are shown in Figure 10. 

Cases 3 and 4 also show similar area distributions and final mass values 
(see Figure 10). And the final structural masses obtained here (15.18 and 
14.94, see Table 10) are lower than the results reported in in Ref. 18 (24.01) 
by more than 35 percent. Between cases 3 and 4 it is seen that case 4 which 
uses 2x2 Riccati equation solutions to generate the startup gains converges 
much faster (see Figure 9). 


8.4 EXAMPLE 3 - ANTENNA STRUCTURE: CDV LINKING (I) 

The third example is an antenna structure consisting of eight aluminium 
beams ( E = 7.3 x 10 6 N /cm 2 , p = 2.77 x 10 3 kg/cm 3 , v = 0.325 ) which have 
thin walled hollow box beam cross sections (see Figure 11). This structure is 
constrained to move vertically (Y - direction) only, so each nodal point has 3 
degrees of freedom (translation, bending and torsion) resulting in the total 18 
degrees of freedom (N = 18). Four translational actuators (A/ = 4) weighing 
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4 kg each are attached to nodes 3, 5, 6 and 7. These actuators are oriented so 
that the force they generate acts in the vertical direction (degrees of freedom 
4, 10, 13 and 16). Two ramp type transient loads are applied to the node 3 
at the same time. One is a vertical force (/, ( 0 ) and the other is a torsional 
moment (f 2 { t)) with respect to finite elements I and 2 which gives anti sym- 
metric exitation. These loads are given as follows: 

f { {()= 333.3 / N, f 2 ( t) = 10.0 x/,( O N • cm 

forO < t < 0.3 seconds, and/, (O=/ 2 (O=0 for t> 0.3 seconds (see Figure 
1 1). Transient responses are considered for the time interval 0 < t < 2 seconds 
and 20 out of 36 complex modes are used to calculate the peak response values. 
The normalizing constants d's tor the right hand eigenvectors (see Eq. (5.8)) 
are chosen to be 10 s . 

Flange and web thicknesses are constrained to be the same, so there are 
three structural design variables for each finite element (B, H and T — T2 — T3, 
see Option 1 1 of Table 1). Structural linking is also used to make the structure 
remain symmetric with respect to the XY plane, which results in the total 15 
independent structural design variables. The initial structure is uniform (B 
= H = 20.0 cm, T = 0.5 cm), and the side constraints are 10.0 cm < B, H 
< 25.0 cm, and 0.1 cm < T < 1.0 cm. 

Move limits of 30 percent for both structural variables and control vari- 
ables are used. Behavior constraints are imposed on: (1) the real part of all the 
retained complex modes ( cr, < -0.5); (2) the fourth and fifth damped fre- 
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quencies (co^ > 8.0 Hz, co ^ > 9.25 Hz ); (3) the peak displacement of nodes 
2. 4, 5 and 7 (| <7, (01 ^ 1.0 cm, i = 1, 7, 10 and 16); and (4) the peak, actuator 
force (| Uj (t) | < 8.5 N,y = 1,2, 3 and 4). 

Initial feedback gains are computed by solving 10 sets of 2 x 2 Riccati 

equations with the control weighting coefficients y's = 1/400 and the state 

weighting matrices [(?,■] = Diag ( c oj, 1 ) i = 1, ... , 10 . In Figure 12, the initial 

closed-loop eigenvalues ( X t = a, + j(a d ) obtained by solving 10 sets of 2 x 2 

Riccati equations are compared with those obtained from a full order Riccati 

equation solution (with weighting matrices [{?] = Diag [ [AT], [A/] ] and 

[/?] = — - — [/] , see Eq. (4.3)), and as can be seen from the plot, these two 
400 

solution methods give almost the same values for the lowest 10 modes. 

Row-wise and column-wise control design variable linking schemes are 
investigated for this example. The maximum number of control design vari- 
ables is relatively large (A/ x 2 N = 144). In this example since M — 4 rather 
than unity, there exist 10 distinct control design variable options. Ten different 
control design variable linking schemes (see Table 2) are imposed on the the 
same initial startup gains from the beginning. 

Iteration histories arc shown in Table 20 and Figure 13 and final struc- 
tural designs are given in Table 21 and Figure 14. Design masses include the 
fixed actuator mass (4x4= 16 kg) as well as the variable structural mass. 

As the freedom in the design space is reduced by imposing more restrictive 
control design variable linking schemes (from case 1 to case 10), it can be 
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clearly seen from the results that: (1) the number of independent control var- 
iables in the optimization loop decreases (from 144 to 1); (2) the optimum mass 
increases (from 163.1 1 kg to 206.06 kg); (3) the total number of iterations de- 
creases and the convergence becomes more robust. 

Even though there is more than 20 percent of difference in the optimum 
mass between case 1 and case 10, all cases show a similar trend in the final 
structural design. Namely, widths and depths of finite elements 1, 2, 5 and 6 
take on their upper bound values and thicknesses of elements 3, 4, 7, and 8 
move to their lower bound values (see Table 21 and Figure 14). 

Final mass values and active constraint sets for the final designs achieved 
in cases 1-10 are summarized in Table 22. In all cases damped frequency, peak 
response and peak control force constraints are active. 

This example problem investigates the effect of row-wise and column-wise 
linking schemes on: (1) the minimum mass achievable; and (2) convergence 
characteristics of the optimization procedure. Generally as more linking is 
imposed, the final mass increases and the design converges faster and more 
smoothly. In Figure 15 optimum masses are compared with the number of 
independent control variables according to the distinct control design variable 
linking options chosen. 
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8.5 EXAMPLE 4 - ANTENNA STRUCTURE: CDV LINKING (II) 


This example is the same as the previous example (Example 3) except that 
here the block type control design variable linking schemes introduced in sec- 
tion 4.4 are investigated. 

8.5.1 Linking on [ H ] 

Cases 1-10 select participation coefficients a 's of [//°p> (see Eq. (4.33)) 
as design variables. These a ' s are further linked so that in each case the 
number of independent control design variables is different. The basic scheme 
employed here is to treat the first (A — 1 ) variables a,, ... , oe A -_ I as independent 
and then link all the remaining variables a K = a K+l = - = a I0 so that the total 
number of independent control design variables after linking is K. For exam- 
ple, when K— 1 there is only one independent design variable after linking, 
since a, = a 2 = — = a )0 . When K = 2 there are two independent design vari- 
ables after linking, namely a, and a 2 = oc 3 — ••• = a I0 . Finally when K~ 10 
there will be ten independent design variables after linking, namely the partic- 
ipation coefficients of the ten basis matrices in Eq. (4.33). 

Iteration histories for cases 1-10 are given numerically in Table 23 and the 
results for cases 1, 2, 4, 6 and 10 are presented graphically in Figure 16. Final 
structural designs are displayed in Table 24 and in Figure 17. 
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8.5.2 Linking on [ Hp ] and [//v] 


Cases 1 1-20 are the same as cases 1-10 except that the position and ve- 
locity parts of the gain matrix are separated. Namely, the ct's of both [//p] (l) 
and [//?]<'■> (see Eq. (4.34)) are candidates for design variables, so that the 
maximum number of independent control design variables after linking is 
doubled from 10 to 20. Iteration histories and final structural designs are 
given in Tables 25-26 and Figures 18-19. 

As the freedom in the design space is increased by choosing more inde- 
pendent control design variables (from case 1 to case 10 and from case 1 1 to 
20), it can be clearly seen from the results that the optimum mass decreases. 
Table 27 summarizes the critical constraint sets ( -0.03 < Gj< 0.0004) at the 
final designs. The fourth damped frequency (cu^) , the peak displacement at 
node 7 (q l6 ) and the peak control force at node 7 (u 4 ) are critical in all 20 cases. 

It is important to note that even with only one or two independent control 
design variables (case 1 and case 1 1), the final objective mass values obtained 
(206.06 kg and 204.16 kg) are about 15% lower than the result reported in 
Ref. 25 (241.97 kg). This can be attributed to the fact that in Ref. 25 the 
control gains are not independent design variables since, for any particular set 
of structural design variables, they are determined from the solution of an 
LQR subproblem. 

In Figure 20 final masses are compared with the number of independent 
control variables (case 1 of Example 3 which has 144 independent CDV s is 
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also shown as a reference solution). Compared to the results for Example 3, 
based on the column-wise or row-wise linking schemes, the new block type 
(basis matrix) linking method (see section 4.4) gives much better results, in the 
sense that with the same or fewer independent control design variables signif- 
icantly lower final design mass can be obtained (compare Figures 15 and 20). 

8.6 EXAMPLE 5 - ANTENNA STRUCTURE: ADDITIONAL 
PROBLEMS 

In this section additional problems using the same antenna structure as 
in Examples 3 and 4 are investigated. All the cases in this section use the same 
startup gains and the same control design variable linking scheme as in case 
10 of Example 4 (namely, startup gains are calculated by solving 10 sets of 2 
x 2 Riccati equations and a,,... , a 10 are independent CDV's, see Eq. (4.33)) 
except for some additional constraints and/or special features. 

8.6.1 Additional Constraints 

Three cases are considered which have additional behavior constraints. 
Case 1 is identical to case 10 of Example 4 except that this case has an addi- 
tional constraint on the control effort such that CE < 20 N 2 • sec. The 4x4 
control weighting matrix [/? CE ] is chosen to be the identity matrix (see Eq. 
(5.21)). Case 2 is identical to case 10 of Example 4 except that in this case 
additional constraints are imposed on the accelerations in the vertical direction 
at nodes 2, 4, 5 and 7 ( | a^(r)| < 100 cml sec 2 , k = 1,7, 10 and 16). In case 
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3 both control effort and acceleration constraints described above are imposed 
at the same time, in addition to all the other constraints. 

Iteration histories and final structural dimensions are given numerically 
in Tables 28-29. and displayed graphically in Figures 21 and 24. Table 36 
shows final mass and critical constraint set results for these three cases. Final 
mass values are increased compared with that for case 10 of Example 4 (170.92 
kg) and this can be attributed primarily to the presence of additional behavior 
constraints. It is interesting to note that in cases 2 and 3, where acceleration 
constraints are considered and found active at the final design, a different type 
of structural design is found, namely one where all the widths (B) take on their 
upper bound values instead of the depths (H) (see Table 29. and Figure 24). 

8.6.2 Updating Fixed Ratios and Truncation of the Gain Matrices 

Cases 4 and 5 are the same as case 10 of Example 4 except that feedback 
gain matrices are modified in order to update the fixed ratios (case 4) or to 
prevent higher uncontrolled modes from being destabilized (case 5). In case 4 
the fixed ratios built into the initial startup gain matrix are updated by solving 
2x2 Riccati equations again at the beginning of each iteration. When 2x2 
Riccati equations are re-solved, weighting matrices are adjusted such that the 
real parts of the closed-loop eigenvalues remain invariant between the updat- 
ing (see Appendix C). In case 5 components of the feedback gain matrix which 
excite uncontrolled higher modes are eliminated by using the method presented 
in Appendix D. 
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Iteration histories and final structural dimensions are given in Tables 
30-31 and Figures 22 and 24. Final masses are similar to each other, differing 
by less than 3 percent difference. In Table 32 the final closed-loop eigenvalues 
arc compared with those of case 10 of Example 4. Most of the uncontrolled 
higher modes (modes 11 to 18) are unstable in case 10 of Example 4, but in 
cases 4 and 5 they are marginally stable. 

8.6.3 Lumped Mass Design Elements 

In this subsection the actuator mass constraints described in section 5.5 
arc included and actuator masses are treated as independent design variables 
(i.e., variable lumped mass design elements). The exponent in Eq. (5.33) is set 
to be unity (c 2 = 1), namely, the relation between the peak control force and 
the required actuator mass is assumed to be proportional with the coefficient 
<V 


Three runs are made with different values of c,. In case 6 the coefficient 
c\ is chosen such that an actuator of 4 kg mass can produce control force of 
8.5 N (i.e., c, = 4 kg j 8.5 N) to provide requirements similar to those in case 
10 of Example 4, where constraints are given on the peak control forces (8.5 
N) for the fixed mass actuators (4 kg). In case 7 more restrictive actuator mass 
constraints are used such that actuators with 8 kg mass are needed to produce 
control forces of 8.5 N (i.e., c, = 8 kg / 8.5 N). And in case 8, c, is chosen such 
that actuators of 2 kg mass are enough to generate control forces of 8.5 N (i.e., 
c x = lkg\ 8.5 N). 
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Iteration histories, final structural and actuator mass values are given in 
Tables 33-35 and Figures 23-24. As expected, case 6 converged to a final de- 
sign similar to that found in case 10 of Example 4, and cases 7 and 8 converged 
to higher and lower final objective mass values, respectively, consistent with 
the higher and lower values used for c { . 

Final mass results for all the cases in this section are summarized along 
with the active constraint sets in Table 36. 


8.7 EXAMPLE 6 - GRILLAGE STRUCTURE: CDV LINKING (II) 

The final example treated here is the 4 by 6 planar grillage structure (see 
Figure 25) which was previously studied in Refs. 9 and 25. It consists of a 
lattice of 10 aluminium frame members placed on 2 foot centers and 
cantilevered trom two fixed supports by 2 foot long flexible beams (E — 
10.5 x 10 6 psi, p = 0.1 lb jin}, v = 0.3). Each solid rectangular member is 2.0 
in wide (fixed) and has an initial depth (variable) of 0.25 in. The members are 
oriented so that the width dimensions lie in the plane of the structure (XZ 
plane). The grillage is modelled using 40 finite elements each of which is 2 foot 
long and the total number of degrees of freedom is 72 (3 per node at 24 nodes). 
A small amount of passive damping (c iVf =0, = 0.00005 , see Eq. (3.2)), 
which gives passive damping ratios between 0.0059 % (1st mode) and 0.36 % 
(20th mode) to the uncontrolled initial structure, is assumed to exist. Four 
torque actuators are placed to provide control torque in the directions as 
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shown in Figure 25. The mass of each actuator (1.296 x 10 -3 lb-sec 2 /'in) is 

modelled as a fixed nonstructural mass. Initial control gains are obtained by 

solving 25 sets of 2 x 2 decoupled Riccati equations with diagonal control 

weighting matrices [/? f ] = — Diag [2, 2, 1, 1] and diagonal 2x2 state 

200 

weighting matrices [£>,] = Diag ( a>}, 1), i = 1 ,... ,25 (see Eq. (4.21)). Initial 
open-loop and closed-loop eigenvalues are given in Table 37. Transient re- 
sponses are considered for a time period 0 < t < 3 seconds and the lowest 40 
out of 144 complex modes are used. 

The total mass is minimized and the grillage is subjected to a transient 
loading at node 3 in the Y direction which is a half sine pulse of magnitude 0.2 
lb and frequency 7 rad/sec. Structural design variable linking is used to im- 
pose symmetry with respect to the XY plane on the structure and as a result 
there are 8 independent structural design variables (see Figure 25). Lower and 
upper limits on the member depths are 0.1 and t.O in, respectively. In this 
example behavior constraints are imposed on: (1) the modal damping factors 
of the first 20 modes (£, > 1%, /= 1,... ,20) ; (2) the transient displacements 
at nodes 1-6, 7, 13, 19, 12, 18, and 24 in the Y - direction (i.e., q{t)< 0.2 in); 
and (3) the transient control torques of all actuators (u(t)< 2.5 lb-in). 

At the initial structural and control design, all the damping ratio and 
control force constraints are satisfied, but transient displacement constraints 
are infeasible by as much as 68 %. Hybrid approximation in terms of depths 
of the members and linear approximation in terms of the control design vari- 
ables are used to generate approximate optimization problems. 
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In this example the number of elements in the feedback gain matrix is 
very large (M x 2/V = 4 x 2-72 = 576), so it is almost impossible to use all the 
gain elements directly as independent design variables. Six cases are investi- 
gated which are similar to cases 1-10 of Example 4, namely, the participation 
coefficients a/s of [//°]<*> (see Eq. (4.33)) are control system design variables. 
These ct's are further linked so that in each case the number of independent 
control design variables is different (see Table 38). 

Iteration histories are given numerically in Table 40 and graphically in 
Figure 26. Final member depths are given in Table 39 and Figure 27. In all 
the cases depth of member 3 (node 13-18) has its lower bound value, and in 
cases 5 and 6 depths of members 5 (node 1-19) and 10 (node 6-24) also have 
lower bound values in addition to member 3. In Figure 28 final mass values 
achieved in each of the six cases are compared with the number of independent 
control design variables used. Similar observations to those made about pre- 
vious examples can be made, namely, as the number of independent control 
design variables is increased from I to 20, the final objective mass value de- 
creases (from 0.1191 to 0.1039, a 12.8 % reduction), but the total number of 
analyses required for the convergence tends to increase. The final closed-loop 
complex eigenvalues and modal damping factors are given in Table 41. 

In all the cases transient displacement constraints at node 1 and the 
transient control force constraints on actuator 3 and 4 are active at the final 
design in addition to the critical damping ratio constraints (noted in Table 41). 
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Chapter IX 

CONCLUSIONS AND RECOMMENDATIONS 
9.1 CONCLUSIONS 

A method has been presented to integrate the design space for 
structural/control system optimization problems in the case of linear state 
feedback control. 

When linear state feedback control is used in the integrated 
structural/control optimization problem, the usual approach is to treat struc- 
tural design variables as the only independent variables in the optimization 
loop. The control gains are subordinated by using the solution to the linear 
optimal control problem to implicitly represent the "optimal" control gains 
corresponding to any given set of structural design variables. In this approach 
the control gains are in effect dependent variables that can be determined for 
any particular structural "plant" once the state and control weighting matrices 
arc specified. In other words, there is no design space freedom on the control 
system design unless the state and control weighting matrices involve some 
candidates for independent design variables which are free to change during 
optimization. While some studies suggest to use coefficients of the weighting 
matrices or control effort as independent design variables in optimization (see 
Refs. 22, 23), this is a rather indirect approach and it does not represent a truly 
integrated formulation of the structural/control system optimization problem. 
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The method presented here treats elements of the feedback gain matrix 
directly as independent design variables together with the structural sizing 
variables and actuator mass variables. In order to reduce the number of in- 
dependent control design variables, the conventional structural design variable 
linking idea is extended to the elements of the feedback gain matrix introduc- 
ing several alternative linking schemes. Also a method for generating effective 
initial startup control gains, which avoids the burden of solving the full order 
Riccati equations (see Section 4.2), is presented. 

The integrated structural/control system design problem is posed as a 
general nonlinear mathematical programming problem. The objective is to 
minimize the total mass of the integrated system. Constraints on dynamic 
stability, damped frequency, control effort, peak transient displacement, ac- 
celeration and control force, and actuator mass are considered. By assuming 
that the external transient disturbances are represented in terms of a truncated 
Fourier series and polynomial terms, all the transient responses and their sen- 
sitivities are derived explicitly in closed form. 

The general optimization problem is solved through the iterative con- 
struction and solution of a sequence of explicit approximate problems based 
on various approximation concepts including new control design variable link- 
ing schemes. Each approximate problem is solved using the feasible direction 
method implemented in CONMIN (Ref. 36) or the modified feasible direction 
method in DOT (Ref. 37). 
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The methodology summarized in this report has been implemented in a 
research computer program, and numerical results for several illustrative ex- 
ample problems have been presented. The number of analyses required to 
achieve convergence (based on a 0.1 percent objective function diminishing 
returns criterion) in most example problems ranged from 8 to 20 and signif- 
icant improvement in the final design objective function values have been 
achieved compared to previously reported results. 

The antenna example was particularly useful as a test-bed for examining 
the alternative control design variable linking schemes introduced here, since 
it involves multiple actuators and moderate model size. These results confirm 
that providing more design freedom by increasing the number of independent 
control design variables (after linking) makes it possible to achieve lower ob- 
jective function values. The trade-off between reduced objective function value 
and increased control variable design freedom is clearly illustrated in Figure 
15 (row-wise and column-wise CDV linking schems) and in Figure 20 (block 
type CDV linking schemes). It should be noted that as the number of inde- 
pendent control variables increases, the number of analyses required for con- 
vergence tends to increase. Although the full extent of the benefit is problem 
dependent, it is often possible to achieve significant reductions in the final ob- 
jective function value by incorporating only a small number of independent 
control design variables into the integrated design space. 

The design variable linking idea has been successfully extended to inte- 
grated structural/control optimization problems, based on full state linear 


feedback, control. In the context of this class of problems, design variable 
linking makes it possible to treat structural design variables (SDV's), actuator 
mass variables and control design variables (CDV's) simultaneously without 
having to deal with prohibitively large numbers of design variables. Further- 
more, the 2 x 2 decoupled Riccati equation solution method makes it possible 
to obtain good startup values for the elements of the gain matrix without the 
computational burden of having to solve a full order (2N x 2 N) Riccati 
equation. The method presented shows promise in the sense that it offers the 
prospect of being able to exploit the benefits of full state feedback and true 
integration without having to: (1) repeatedly solve large 2 N x 2 N Riccati 
equations; (2) deal with extremely large numbers of independent design vari- 
ables. 

9.2 RECOMMENDATIONS 

The method presented in this work can be used to optimize frame/truss 
structures augmented by a linear state feedback controller. While this work 
represents an important step towards the goal of integrating the structural and 
the control system design processes, several improvements and extensions can 
be made which will lead to increased efficiency and broader applicability. 

Since the majority of the time for each design iteration lies in the analysis 
(analysis and sensitivity analysis), reducing the total number of analyses re- 
quired for the convergence will lead to significant reductions in the total cost. 
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The key to fewer design iterations is the construction of more accurate ap- 
proximate optimization problems so that larger move limits on the design var- 
iables can be used during each stage. One way to improve the accuracy of the 
approximations is to use intermediate response quantities such as modal ener- 
gies when constructing explicit approximations for the real (a) and imaginary 
{<x> d ) parts of complex eigenvalues (see Ref. 9). 

Another important factor in the total cost of optimization is the time re- 
quired for each analysis. If the order of the model is reduced by using appro- 
priate basis vectors (for example natural modes of the structure), the analysis 
cost can be significantly reduced. 

Currently a full order finite element model of the structure is used for the 
controller design (in other words, elements of the feedback gain matrix, [//] , 
are optimized in the physical coordinate system) assuming all the dynamic 
displacements and velocities are directly available. However, for practical 
purposes an observer or a state estimator is required to reconstruct the state 
values from the sensor outputs. When the observer is constructed for the full 
order model the task is rather difficult, since it can increase the order of the 
entire set of dynamic equations up to AN. Once the original model is reduced 
as mentioned above, building the observer becomes much more tractable. 
When an observer is used, the same linking idea can be applied to the observer 
gains. It should be noted that when the observer gains are also used as inde- 
pendent design variables, various linking schemes can be used for the observer 
gains. Furthermore, it will be necessary to consider additional constraints, for 
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example the real parts of the observer poles should have values that are more 
negative than the real parts of the controller poles (this will make the estimator 
errors decay faster than the controlled states). 

In Appendix C a method is presented for updating the fixed ratios be- 
tween gain matrix elements (associated with a particular linking scheme), while 
constraining the real parts of the closed-loop eigenvalues (a ; s) to be unchanged 
during the update. When modal damping ratios (c* 5 ) are constrained rather 
than o’s , the possibility of a similar updating scheme needs to be investigated. 

Other areas for further investigation which will be important for the sol- 
ution of more realistic problems include consideration of: (1) stress constraints 
(both static and dynamic); (2) multiple loading conditions; (3) overdamped 
modes (pure real closed-loop eigenvalues); (4) sensor/actuator location issues 
and/or observability and controllability requirements; and (5) sensor/actuator 
time delay effects. 
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Appendix A 


ELEMENT STIFFNESS AND MASS MATRICES FOR SPACE FRAME 

ELEMENT 

For a 3 dimensional frame element shown in Figure 1 the vector of the 
e-th clement nodal degrees of freedom in the local coordinates is 

{q e } T = [«!, v,, w,, e xl , e y[ , 9 zl , u 2 , v 2 , W 2 , e x2 , e y2 , Q zl } (-4.1) 

where «, v, and w are displacements along x, y, and z axes, d x , d y and d 2 are 
rotations about x, y, and z axes, respectively, and subscripts 1 and 2 represent 
each end point 1 and 2. Then the element stiffness matrix is 

AO 0 000 -AO 0 0 00 

1 21 J l 2 0 0 0 61 J l 0 - 12 Iji 1 0 0 0 61 J l 

12 y/ 2 o -6 y/ ooo - 12 y/ 2 0 -6 y/ 0 

GJjE 0 0 0 0 0 - GJ\E 0 0 

4/j, 0 0 0 6 y/ 0 2I y 0 

4/y 0 6 I y jl 0 0 0 2I Z 

A 0 0 0 0 0 

12/,// 2 0 0 0 - 61 J l 

Sym. 12 y/ 2 0 6 l y H 0 

GJjE 0 0 

4/y 0 

4 /, 
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where E is the Young's modulus, G is the shear modulus, A is the area, / is the 
length, GJ is the torsional rigidity, and I y and I 2 are the area moment of inertia 
of the cross section with respect to the centroidal principal axes y and z, re- 
spectively. 

The element mass matrix is 

140 00 0 00 70 00 0 00 

156 0 0 0 22 / 0 54 0 0 0 13 / 

156 0 - 22 / 0 0 0 54 0 13/0 

140 /p/ /I 0 0 0 0 0 10I p IA 0 0 

4/ 2 0 0 0 - 13 / 0 - 3/ 2 0 

4/ 2 0 13 / 0 0 0 - 3/ 2 

140 0 0 0 0 0 

156 0 0 0 22 / 

Sym. 156 0 22 / 0 

140 1 p l A 0 0 

4/ 2 0 

4/ 2 

where p is the mass density and I p is the polar area moment of inertia of the 
cross section. 
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Appendix B 

SOLUTION OF 2 x 2 RICCATI EQUATION 


Equation (4.23) can be written in equivalent form as 


P\ 1 

5 1 2 

5i 2 

5 22 

5i 1 

5 12 

5l2 

5 22 


0 1 


0 — coj 


P\\ 

P 12 


fill 

0 

2 

<°i - c i\ 

+ 

1 - c i 


Pi 2 

P 22 

+ 

0 

Q 22 


(5.1) 


[ 0 ] 





512 
5 22 


= [0] 


where the 2 x 2 matrix [5 ( ] is positive definite (i.e., p‘ u > 0, p‘ 22 >0 and 
P\ i Pn ~ iPnf > 0) and Q l u , Q‘ n > 0, y t > 0. 


Since the above equation is symmetric, there are 3 independent scalar 


equations as follows: 

Wi(p[ 2 ) 2 + 2 c o}p\ 2 - eij = o (B.2.a) 

W'iPn) 2 + 2 Cj p '22 - G 22 - 2 P\2 = 0 (B.2.b) 

^ (512522 — P\\ + c iP\2 + P22 ~ ® (5.2. c) 

where 

If; = }|Vi) r MW r {»i) > 0 (S.3) 

' ( 
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From Eq. (B.2), 


, _ -« ? ± VW + wygn 

P12 ■ HA 


7 ^+ ||/ , ej 2 + 2 »>| 2 

F22 - 


(5.4. a) 


(5.4.6) 



Since />$ 2 > 0 and q > 0 the minus sign in front of the square root in Eq. 
(B.4.b) is dropped and since p' u > 0 the minus signs of Eq. (B.4.c) disappear 
so that a unique solution for [5,] can be determined as follows: 


512 


— (vj + \ ; cu^ + 


(B.5.a) 



(5.5.6) 
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Appendix C 

UPDATING FIXED RATIOS OF CONTROL GAINS 

As mentioned in Section 4.2, solving for the initial feedback gain matrix 
and imposing some kind of linking on the feedback gains fixes the relative va- 
lues of those elements through the final stage. When the decoupled Riccati 
equation solution method presented in subsection 4.2.2 is used, this problem 
can be relaxed. 

When the 2x2 Riccati equations are solved again for a different struc- 
ture, there should be some rule to choose weighting matrices [(7J and [/?,] of 
Eq. (4.21). Considering that the real parts of the eigenvalues ( o's ) play an 
important role on the dynamic responses as well as themselves are the behavior 
constraints, the scheme presented here has focused on the continuity of the 
o's during the resolution of the feedback gain matrix. 

First, the further assumption is made that the coupling effects of Eq. 
(4.27) on the closed-loop eigenvalues are negligible. Then neglecting the right 
hand side of Eq. (4.27), the i-th equation for the complex eigenvalue X t becomes 

X 2 t + (c,- -I- W t P 22 ) + WiP\l) = 0 (£-1) 

where 

w i = T- [*] r { V, } (C.2) 

' / 
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Then the i-th closed-loop eigenvalue X- t becomes (assuming underdamped mo- 
tion) 


x i = ± 

- (c f + W, P 22 ) ± 7\/ 4 ((U^ + ( c i + ^'^ 2 2 ) 2 

= 2 

From Eq. (C.2)-(C.5) noting that { v ( } , [6], cu? and c, are fixed for a given 
structure, the complex eigenvalues X t can be assigned arbitrarily by adjusting 
IF/ , p\ 2 and p‘ 22 or equivalently y ; , <2i 1 and Q 22 . 


The ratio between the state weighting matrix [(>,] = Diag (£>' u , Q[ 2 ) and 
the control weighting matrix [/?,] = y, [/] determines the relative magnitude 
of states and control forces. This means that the ratio between [(?,] and [/?,] 
(or Q\ u Q { 2 and y f ) will determine the damping effect in the i-th mode, so by 
changing only y, the real part of the closed-loop eigenvalue can be assigned 
(within some bounds). 


The foregoing observations suggest that for some iterations (for example 
for every K iterations) the fixed ratios within the feedback gain matrix can be 
updated by changing the y-s in a manner which forces the real part of the 
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closed-loop eigenvalues to coincide with the approximate real part of the ei- 
genvalues ( o' s ) from the previous iteration, i.e. from Eqs. (C.5) and (C.4) 

c i + W iP22 

a = _ 

' 2 

^ == _ (C.6) 

= - yV f + "'iQ'n - 2coj + 2 Jw* + 

or 

w , = yy-y (* + 2 ej, - 2 Jq‘„x + (&,) 2 + (d 2 ) 2 ) (C. 7) 

\^22/ 

where 

-V = <>j 2 (4 ? • + 2cu- - c, 2 ) 

Then from Eq. (C.2) 

?, = 1 y(v,) 7 'OTM 7 'f» i } (C.8) 

With a new set of y/s the decoupled Riccati equations are solved again so 
that the relative values of the elements in the gain matrix are updated accord- 
ing to the changing structure preserving continuity of the real parts of the 
complex eigenvalues between the iterations. When all the y's converge to the 
previous values, this updating option can be turned off. 
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Appendix D 

ABOUT SPILLOVER EFFECT 


In this appendix the existence of spillover effect (see Ref. 38) of the for- 
mulation used in this study is investigated. Consider Eqs. (4.1) and (3.9) again 
for convenience 

[AO { ? 1 + [C] { q } + [fl { q > = [4] { u } (0.1) 

{ u } =-[//,]{?}- [//„] { q ) (0.2) 

The nodal degree of freedom rector, { q } , can be written without any trun- 
cation as follows (compare with Eq. (4.9)): 

{q} = ( z c) + C v ul{*u) 

where { z c } is the r x 1 controlled normal coordinate vector, [F c ] is the N x r 
eigenmatrix consisting of the controlled r normal eigenvectors, { z u } is the 
( V - r) x 1 uncontrolled normal coordinate vector, and [VJ] is the 
N x (A' — r) eigenmatrix corresponding to uncontrolled modes { z v } . 

Eigenmatrices [F c ] and [VJ\ in Eq. (D.3) can be normalized such that 

[t'cfw ] [f'ci = m, [i / u ] 7 [.v/] cm = m 

(D.4) 

[> / d r [AO mt'] = CO], [K/[M] [f'd = [0] 

and 
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[^c/rc] [K c ] = [C c ], [t' c ] r [K] [T c ] = [A c ] 

[ (/ t'] r [C][^ l ,] =[C y ], [l / y] r ra[F t ,] =[Ay] <Z) ' 5) 

Substituting Eq. (D.3) into Eq. (D.l) and premultiplying by [F c ] r and 
results in 

^c) + [Cd(^c} + [Ad { z c } = lV c f[bl{u} (D. 6) 

+ C c d ( *u ) + [ A d ( z u ) = [^d r [d{«} {D.l) 

On the right hand side of Eq. (D.7), [Fd r | >] # [0], which means there ex- 
ists control spillover. 

Substituting Eq. (D.3) into Eq. (D.2) leads to 

! " ) = - (t*0J ( z C ) + W U'd ! i c i) 

- flyy c y v i ( ) + w cm f >) (Z> ' 8) 

On the right hand side of Eq. (D.8), those matrices in front of the uncontrolled 
modes are not zero (i.e., [//,] [FJ # [0], [// v ] [FJ * [0] ) which can be 
interpreted as observation spillover. 

As can be seen from Eqs. (D.7) and (D.8), there exist both control and 
observation spillover. According to Ref. 38, the closed-loop system has po- 
tential instability when the system has both observation and control spillover 
and it is more important to eliminate observation spillover in order to avoid 
instability. Therefore, the desired control input, { u * } , should contain only 
controlled modes, { z c } , namely 
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{«*) =- \.n p ] [v c l f z c! - Weenie) 

Premultiplying [P C ] 7 '[A/] to Eq. (D.3) and using relations in Eq. (D.4) results 


in 


{z c } = IVcfim {q} 


(DAO) 


Substituting Eq. (D.10) into Eq. (D.9) gives 

{u*} = - [Hp] [ Vq\ [ [A/] { q } - 

= - VHp] { q } - [//*] { <7 } 


where 

[«;] = IH^lVcUVdfW 

and 

[<] = 


(ZD. 1 2) 


(D. 13 ) 


In summary, destabilization of uncontrolled higher modes can be pre- 
vented by using the truncated feedback gain matrices shown above since these 
matrices are orthogonal to the uncontrolled higher modes [P^] and will elimi- 
nate observation spillover. 
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Appendix E 

CLOSED FORM SOLUTION OF TRANSIENT ANALYSIS 
E.l CLOSED FORM RESPONSE CALCULATION 


Eq. (5.17) can be written in explicit closed form as 


N, 


a 


/V p 


'/,(') = { >K) T IE] I £ ( IC k{FC} k + IS k { FS ) k ) + £ fP p {FP} p 

\* =1 p= 0 J 

+ e ' 1 1 7, (0) when 0 ' (^-0 

*7/(0 = when t>t 


7 


7 


where 


/C t 




cosQ* r <7 t 


2 . ,2 

i 


+ a 


+ O* sin ^ cosQ*r) 


AS* = f e A ' (, - T) sinQ*T^T 

Jn 


2 . ,2 

i 


n* + a 


(f2*e Ai/ - O* cosQ*r — A,- sinO*/) 


IP = [ e X ^~ z) x p dx 

Jn 


(£. 2 ) 


(£. 3 ) 


(£. 4 ) 


When p — 0 (which corresponds to a step input), 
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(£. 5 ) 


IP Q = [ e A,( ' J) dx = - l) 

4 A /' 


and when p = 1 (which corresponds to a ramp input), 


ip i 3 f e ' 1 ' 


(f-i) , 1/1 A, / , 1 

trfr = Tlr 


(£. 6 ) 


The above expressions can be derived by integration by parts. It should be 
noted that for a given set of design variables, r\ t {t f ) needs to be calculated 
only once and for t > t f the second formula of Eq. (E.l) is used. 


E.2 CLOSED FORM RESPONSE SENSITIVITY CALCULATION 


Differentiating Eq. (E.l) with respect to a yields 


Wi(Q 


d{*h} 


T 


da 


da 


[£] + {' M 


T 3[£] 
da 


N 


N P 


x 


V ( !C k {FC) k + is k {FS} k ) + Y, IPpiPP}, 

P= 0 


*=1 


+ {*/}[£! 


k= 1 V ' P= 0 


(E.l) 


dk: 1 , 5 ^( 0 ) 

+ + * ' —i — . 

da da 


when 0 < t < t 


7 


= — (r- 7e A,( '“ 7> ^ ; M + e A <'7 7 M ll when 

da da ; ; dot 
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where 


cICi 


dOL 


= f e^ 1 ^cosQ^r dx 

SO. Jg * 

= ^L[ 

da Jr 


dXi Ct e M^) (f _ 
0 


dX, 


da Cll + xj 


(t - t) cos Cl/c x dx 

(t Xj /' ' + e k ‘ { - cos Q k t - 2X f /C*) 


cIS L 


da 


J e^‘ T ^sin D. k xdx 

17 sin£ V* 

«, i 


5a + xj 


(t n k e X ‘ 1 - sin Q^r - 2^/S*) 


5//>n 


5a 


5/P, 


(?a 


dot J 

e 

0 

dX { 

da , 

r* 

J o 

dXj 

( 

da ' 

r 

J_| 

V , 
/ 

& 

5a J 

c 

0 

dX; 

c* 

da . 

'0 

5A ; - 

1 


t(t-r) 


dx 


Xit 


■i(t-x) 


t 


X; 


x dx 


X, 


d(X xj 


l L 


i x dx 


' + rir' + i 


Again the above expressions can be derived using integration by parts. 


(£. 8 ) 


(£.9) 


(£. 10 ) 


(£. 11 ) 
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TABLE 1 


Choice of Structural Design Variables 


Option 

Design 

Variables 

Linking within 
an Element 

i 

B and II 

NA 

2 

II 

B =dl 

3 

B 

NA 

4 

11 

NA 

5 

B, II, T2 and T3 

NA 

6 

T2 and T3 

NA 

7 

T2 and T3 

B = a T3, I I = 0 T2 

8 

T3 

T2 = a T3 

9 

T2 

NA 

10 

T3 

NA 

11 

B, II and T3 

T2 = a T3 

12 

II, T2 and T3 

B = a H 

13 

II and T3 

B = a H, T2 = 0 T3 


a, /? : preassigned constants 

NA : no linking within a cross section. 

n.b. : options 5-13 apply to the box beam type element only, (see Figure 1) 
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TABLE 2 


Control Design Variable Linking Options 


Option 

Description 

Design 

Variables 

No. of 
D.V.'s 

i 

totally unlinked 

elements of [//] 

M x 2N 

2 

[//,] fixed, 
[//„] unlinked 

elements of [//J 

M x N 

3 

columns of [//] 
linked 

coefficients of 
columns of [//] 

2N 

4 

C/43 fixed, 

columns of [//,] linked 

coefficients of 
columns of [//„] 

N 

5 

rows of [//,] 
and [//„] linked 

coefficients of 
rows of [//J and [//J 

2M 

6 

rows of [//] 
linked 

coefficients of 
rows of [//] 

M 

7 

[f/] fixed, 
rows of [//,] linked 

coefficients of 
rows of [//„] 

M 

8 

[//,]. [//,] 
linked 

coefficients of 

ra. m 

2 

9 

[f/J fixed, 
[//,] linked 

coefficient of [//J 

1 

10 

[//] linked 

coefficient of [//] 

1 


[//] : M x 2N feedback gain matrix 

[//] : position part of [//] ( M x N ) 

C/4] : velocity part of [//] ( M x N ) 
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TABLE 3 


Iteration History for Example I: Cantilever Beam, Cases 1-5 
Control Design Variable Linking (I) 


Analysis 

Case 1 

i 

2141.6 

2 

1351.1 

3 

833.30 

4 

595.18 

5 

502.20 

6 

461.45 

7 

444.36 

8 

428.04 

9 

427.74 

10 

423.71 

11 

423.70 

12 

423.67 

13 


14 


15 


16 



Total Mass(kg) 


Case 2 

Case 3 

2141.6 

2141.6 

1352.1 

1351.4 

841.33 

841.18 

896.81 

598.67 

621.41 

502.45 

519.07 

460.13 

479.59 

435.03 

467.22 

428.93 

463.20 

428.93 

455.93 

428.93 

454.28 


454.18 


453.07 


452.20 


451.95 


451.94 



Case 4 

Case 5 

2141.6 

2141.6 

1453.5 

1453.6 

956.15 

956.45 

658.18 

658.57 

533.76 

534.30 

485.52 

486.13 

472.54 

473.19 

472.54 

473.19 

472.54 

473.19 


TABLE 4 

Final Cross Sectional Dimensions for Example 1: Cantilever Beam, Cases 1-5 
Control Design Variable Linking (I) 


Final Cross Sectional Dimensions(cm) 



Case 1 

Case 2 

Case 3 

Case 4 

Case 5 

12 

T3 

m 

0.6356 

0.5000° 

0.5290 

0.5000° 

0.7310 

0.5000“ 

0.7340 

0.5000° 


9 indicates lower bound value 
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TABLE 5 


Summary of Example 1: Cantilever Beam, Cases 1-5 
Control Design Variable Linking (I) 



Case 1 

Case 2 

Case 3 

Case 4 

Case 5 

CDV 

Linking 

Option* 

1 or 3 

2 or 4 

5 or 8 

7 or 9 

6 or 10 

No. of 
Total 
Analyses 

12 

16 

10 

9 

9 

Final 

Active 

Constraints 

co,, <7 w 

(JL> d U 

(o d q u 

o’ io 

O,o 

Number 

of 

CDV 

40 

20 

2 

1 

1 

Final 

Mass 

(kg) 

423.67 

451.94 

428.93 

472.54 

473.19 


+ : See Table 2 
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TABLE 6 


Analysis 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


History for 

Example 1: Cantilever Beam, Cases 6-9 

Different Initial Feedback 

Gains 



Total Mass(kg) 



Case 6 

Case 7 

Case 8 

Case 9 

2141.6 

2141.6 

2141.6 

2141.6 

1891.0 

1433.9 

2151.5 

1559.4 

927.12 

1035.4 

1644.0 

1091.8 

549.12 

774.95 

1 170.6 

780.60 

512.22 

595.87 

836.46 

582.35 

465.22 

491.20 

622.14 

456.71 

437.93 

458.57 

489.56 

435.84 

429.40 

431.84 

444.44 

428.09 

429.40 

431.84 

434.64 

427.81 

429.40 

431.84 

434.64 

428.63 



434.64 

428.63 




428.63 


TABLE 7 

Final Cross Sectional Dimensions for Example 1: Cantilever Beam, Cases 6-9 


Different Initial Feedback Gains 

Final Cross Sectional Dimensions(cm) 



Case 6 Case 7 Case 8 Case 9 

T2 

T3 

0.5311 0.5424 0.5554 0.5276 

0.5000° 0.5000° 0.5000“ 0.5000° 


° indicates lower bound value 
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TABLE 8 


Summary of Example 1: Cantilever Beam, Cases 6-9 
Different Initial Feedback Gains 



Case 6 

Case 7 

Case 8 

Case 9 

Initial Gains 

m = [0] 

direct output 

Full Riccati 

2x2 Riccati 

No. of 
Unlinked 
Iterations 

3 

3 

3 

3 

No. of 
Total 
Analyses 

10 

10 

11 

12 

Final 

Active 

Constraints 

w d q u 

w d q u 

a w (o d q u 

w d q u 

Final 

Mass 

(kg) 

429.40 

431.84 

434.64 

428.63 
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TABLE 9 


Nodal Point Coordinates for Example 2: ACOSS FOUR 


Node 

X 

Y 

Z* 

i 

0.0 

0.0 

8.165 

2 

-5.0 

-2.8868 

0.0 

3 

5.0 

-2.8868 

0.0 

4 

0.0 

5.7735 

0.0 

5 

-6.0 

-1.1547 

-2.0 

6 

-4.0 

-4.6188 

-2.0 

7 

4.0 

-4.6188 

-2.0 

8 

6.0 

-1.1547 

-2.0 

9 

-2.0 

5.7735 

-2.0 

10 

2.0 

5.7735 

-2.0 


* : the origin (0,0,0) lies in the plane of nodes 2,3 and 4 
(2 units shift in Z coordinate origin from Ref. 1 1) 


104 





TABLE 10 


Initial and Final Cross Sectional Areas of Truss Elements for Example 2: 

ACOSS FOUR 


Element 

(Nodes) 

Initial 
A and B 

Final 
Case 1 

Final 
Case 2 

Final 
Case 3 

Final 
Case 4 

1 (1-2) 

1000. 

258.2 

251.7 

258.6 

290.4 

2 (2-3) 

1000. 

164.8 

162.6 

260.1 

243.8 

3(1-3) 

100. 

153.3 

149.3 

168.5 

134.9 

4(1-4) 

100. 

152.4 

149.4 

101.4 

131.8 

5 (2-4) 

1000. 

165.8 

162.6 

243.5 

245.0 

6 (3-4) 

1000. 

282.8 

274.5 

291.2 

286.9 

7 (2-5) 

100. 

162.4 

156.4 

76.2 

78.9 

8 (2-6) 

100. 

161.6 

156.5 

137.9 

72.8 

9 (3-7) 

100. 

102.9 

82.0 

110.3 

98.4 

10 (3-8) 

100. 

219.8 

225.1 

102.1 

110.2 

1 1 (4-9) 

100. 

103.8 

81.9 

94.2 

99.6 

12 (4-10) 

100. 

219.6 

225.1 

166.8 

107.9 

Structural 

Mass 

43.697 

14.518 

14.124 

15.177 

14.935 

Control 

Design 

Variable 

1.0 

5.2283* 

2.8781** 

2.9479* 

1.5385** 

Number of 
Analyses 
Required 

- 

11 

12 

13 

8 


* : represents the ratio of the final feedback gain with respect to 
the initial gain of Design A (See Table 1 1) 

+ * : represents the ratio of the final feedback gain with respect to 
the initial gain of Design B (See Table 13) 
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TABLE 11 


Feedback Gains of Initial Design A for Example 2: ACOSS FOUR 
(Feedback Gain Calculated by Solving Full Order Riccati Equations 
with Weighting Matrices [ Q~\ = \_f] 2N and [/?] = 


Position Gains 


Node Direction 

Actuator 1 

Actuator 2 

Actuator 3 

Actuator 4 

Actuator 5 

Actuator 6 


X 

-0.17112 

-0.15531 

0.04389 

-0.00572 

0.01124 

-0.02256 

i 

Y 

-0.08062 

-0.10785 

-0.01233 

-0.02279 

0.04415 

0.00646 


Z 

-0.23050 

-0.23037 

-0.00667 

-0.15163 

-0.00672 

-0.15156 


X 

0.50329 

0.35341 

-0.26035 

0.19237 

-0.15928 

0.15346 

2 

Y 

0.11764 

0.37713 

-0.03369 

0.06621 

-0.20857 

0.13342 


z ! 

0.27425 

0.27412 

-0.14808 

0.04077 

-0.14804 

0.04074 


X 

-0.37381 

-0.12643 

0.07734 

-0.17505 

-0.02042 

-0.01554 

3 

Y 

0.21863 

0.23719 

0.14667 

-0.09655 

0.07495 

-0.13230 


Z 

0.07594 

0.03240 

0.08307 

0.05708 

0.04116 

-0.02841 


X 

0.14220 

0.00236 

0.05463 

-0.12230 

0.16562 

-0.17110 

4 

Y 

-0.22824 

-0.43309 

-0.05519 

0.05272 

-0.00637 

-0.10331 


Z 

0.03239 

0.07595 

0.04116 

-0.02840 

0.08307 

0.05707 

Velocity Gains 


X 

0.33432 

-0.01836 

-0.23952 

-0.53569 

0.34669 

-0.30211 

1 

Y 

-0.21421 

0.39664 

0.53864 

-0.03958 

-0.47675 

-0.44411 


Z 

-0.19287 

-0.19284 

0.16896 

0.08101 

0.16895 

0.08103 


X 

0.84740 

-0.10095 

0.01224 

0.01220 

0.12992 

-0.09151 

2 

Y 

-0.60568 

1.03660 

0.14295 

-0.11270 

-0.06087 

0.06691 


Z 

1.12920 

1.12920 

-0.01356 

0.01441 

-0.01355 

0.01440 


X 

-0.08273 

-0.16078 

0.47486 

-0.52061 

0.17640 

-0.17019 

3 

Y 

-0.09619 

0.20450 

0.85744 

-0.81356 

0.05456 

-0.05683 


Z 

-0.00656 

0.00742 

0.84300 

0.82788 

0.01139 

0.01026 


X 

0.09672 

-0.12468 

0.13546 

-0.13433 

0.97992 

-0.96481 

4 

Y 

-0.24152 

-0.02350 

0.12552 

-0.11899 

-0.01753 

-0.04402 


Z 

0.00741 

-0.00656 

0.01139 

0.01026 

0.84299 

0.82788 
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TABLE 12 


Natural Frequencies, Closed-Loop Eigenvalues and Modal Damping Ratios 
for Initial Design A in Example 2: ACOSS FOUR 


Mode 

Number 

Natural 

Frequency 

(<u 2 ) 

Closed-loop 
Eigenvalue 
(<T ±J(O d ) 

Damping 

Ratio 

to 

1 

1.8010 

-0.05196 + j 1.34173 

0.03869 

2 

2.7715 

-0.07716 ± j 1.66395 

0.04632 

3 

8.3563 

-0.15106 ± j 2.88766 

0.05224 

4 

8.7465 

-0.16782 ± j 2.95366 

0.05673 

5 

11.548 

-0.20178 ± j 3.39314 

0.05936 

6 

17.678 

-0.25698 ± j 4.19752 

0.06111 

7 

21.735 

-0.25094 ± j 4.65595 

0.05382 

8 

22.613 

-0.24350 ± j 4.74958 

0.05120 

9 

72.923 

-0.20643 ± j 8.53709 

0.02417 

10 

85.574 

-0.19523 ± j 9.24860 

0.02110 

11 

105.78 

-0.15129 + j 10.28381 

0.01471 

12 

166.55 

-0.05871 + j 12.90510 

0.00455 
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TABLE 13 


Feedback Gains of Initial Design B for Example 2: ACOSS FOUR 
(Feedback Gain Calculated by Solving 12 Sets of 2 x 2 Riccati Equations 
with Weighting Matrices [<2J = Diag (c oj, 1) and [/?,] = [T ] M ) 


Position Gains 


Node Direction 

Actuator 1 

Actuator 2 

Actuator 3 

Actuator 4 

Actuator 5 

Actuator 6 


X ! 

-0.00027 

-0.00040 

-0.00009 

-0.00015 

0.00008 

-0.00011 

1 

Y 

-0.00031 

-0.00008 

0.00015 

-0.00004 

-0.00015 

- 0.0001 1 

Z 

-0.00078 

-0.00078 

-0.00005 

0.00001 

-0.00005 

0.00001 


X 

0.17636 

-0.17577 

-0.00057 

0.00059 

-0.00008 

0.00009 

2 

Y 

-0.30476 

0.30509 

0.00024 

-0.00023 

-0.00062 

0.00063 


Z 

0.35224 

0.35224 

0.00014 

-0.00019 

0.00014 

-0.00019 


X 

-0.00044 

0.00087 

0.17609 

-0.17614 

-0.00008 

0.00013 

3 

Y 

-0.00013 

0.00026 

0.30499 

-0.30492 

0.00070 

-0.00069 


Z 

-0.00003 

-0.00003 

0.35173 

0.35169 

0.00010 

0.00006 


X 

0.00066 

-0.00033 

0.00057 

-0.00053 

0.35215 

-0.35211 

4 

Y 

0.00062 

-0.00031 

-0.00042 

0.00046 

- 0.0000 1 

-0.00007 


Z 

-0.00003 

-0.00003 

0.00010 

0.00006 

0.35173 

0.35169 

Velocity Gains 


X 

0.48833 

-0.05431 

-0.35118 

-0.87244 

0.55421 

-0.50961 

i 

Y 

-0.34465 

0.59525 

0.84274 

-0.08477 

-0.72551 

-0.71314 

Z 

-0.39794 

-0.39787 

0.30804 

0.11272 

0.30801 

0.11277 


X 

1.63810 

-0.19273 

0.00299 

0.07450 

0.21374 

-0.13354 

2 

Y 

-1.16800 

2.00250 

0.24509 

-0.19719 

-0.11992 

0.16311 


Z 

2.18680 

2.18680 

-0.03887 

0.05238 

-0.03886 

0.05237 


X 

-0.22706 

-0.30616 

0.94046 

-0.99591 

0.30379 

-0.29970 

3 

Y 

-0.16336 

0.36588 

1.64910 

-1.58520 

0.14310 

-0.11673 


Z 

0.00539 

0.00812 

1.62040 

1.59280 

0.03972 

0.01485 


X 

0.16378 

-0.25501 

0.27585 

-0.25098 

1.89820 

-1.87060 

4 

Y 

-0.44813 

-0.11487 

0.19162 

-0.20122 

-0.01018 

-0.06976 


Z 

0.00811 

0.00541 

0.03972 

0.01484 

1.62040 

1.59280 
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TABLE 14 


Natural Frequencies, Closed-Loop Eigenvalues and. Modal Damping Ratios 
for Initial Design B in Example 2: ACOSS FOUR 


Mode 

Number 

Natural 

Frequency 

(co 2 ) 

Closed-loop 
Eigenvalue 
(<f ±ju>d) 

Damping 

Ratio 

(O 

i 

1.8010 

-0.08264 ± j 1.34786 

0.06120 

2 

2.7715 

-0.13181 ± j 1.67555 

0.07843 

3 

8.3563 

-0.28566 ± j 2.90017 

0.09802 

4 

8.7465 

-0.31846 ± j 2.96470 

0.10680 

5 

11.548 

-0.38801 ± j 3.40683 

0.11316 

6 

17.678 

-0.49915 ± j 4.20444 

0.11789 

7 

21.735 

-0.49119 ± j 4.64576 

0.10514 

8 

22.613 

-0.47814 ± j 4.72990 

0.10058 

9 

72.923 

-0.41002 ± j 8.50867 

0.04813 

10 

85.574 

-0.38734 ± j 9.21960 

0.04198 

11 

105.78 

-0.30212 ± j 10.26629 

0.02942 

12 

166.55 

-0.11446 + j 12.88505 

0.00888 
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TABLE 15 


Natural Frequencies, Closed-Loop Eigenvalues and Modal Damping Ratios 
for Final Design of Case 1 in Example 2: ACOSS FOUR 


Mode 

Number 

Natural 

Frequency 

K) 

Closed-loop 
Eigenvalue 
0 ±jw d ) 

Damping 

Ratio 

(0 

i 

1.7917 

-0.22511 ± j 1.34857* 

0.16465 

2 

2.5221 

-0.16755 + j 1.59354+ 

0.10456 

3 

7.7894 

-0.45422 ± j 2.90936 

0.15426 

4 

8.3369 

-0.171 10 ± j 2.91693 

0.05856 

5 

13.329 

-0.26529 ± j 3.72512 

0.07104 

6 

21.781 

-1.31127 ± j 4.39815 

0.28571 

7 

24.304 

- 1.06505 ± j 4.73268 

0.21955 

8 

33.428 

- 1.64435 ± j 5.54856 

0.28414 

9 

37.101 

-0.32458 ± j 6.05921 

0.05349 

10 

44.278 

— 1.72945 ± j 6.22344 

0.26775 

11 

45.840 

-1.47318 + j 6.40279 

0.22423 

12 

48.746 

— 1.69234 ± j 6.67545 

0.24574 


+ : critical constraints 
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TABLE 16 


Natural Frequencies, Closed-Loop Eigenvalues and Modal Damping Ratios 
for Final Design of Case 2 in Example 2: ACOSS FOUR 


Mode 

Number 

Natural 

Frequency 

(* 2 ) 

Closed-loop 
Fi gen value 
(<T ±j(O d ) 

Damping 

Ratio 

(0 

1 

1.6850 

-0.22690 + j 1 .3343 1 * 

0.16764 

2 

2.4551 

-0.15342 ± j 1.59436* 

0.09579 


6.8799 

-0.62833 + j 2.90255 

0.21157 

4 

8.1448 

-0.16614 ± j 2.92190 

0.05677 

5 

13.008 

-0.25024 ± j 3.71749 

0.06716 

6 

18.175 

-1.38884 ± j 3.96390 

0.33066 

7 

22.033 

-1.09006 ± j 4.23768 

0.24912 

8 

32.395 

-1.76219 ± j 5.42404 

0.30899 

9 

35.531 

-0.21716 ± j 5.96693 

0.03637 

10 

44.649 

-1.92262 ± j 6.05257 

0.30275 

1 1 

44.946 

-1.57552 ± j 6.45964 

0.23696 

12 

49.007 

-1.81167 ± j 6.66486 

0.26231 


: critical constraints 





TABLE 17 


Natural Frequencies, Closed-Loop Eigenvalues and Modal Damping Ratios 
for Final Design of Case 3 in Example 2: ACOSS FOUR 


Mode 

Number 

Natural 

Frequency 

K) 

Closed-loop 
Eigenvalue 
(cr ±ju d ) 

Damping 

Ratio 

(O 

i 

1.7695 

-0.14827 ± j 1.33906* 

0.11006* 

■> 

2.2701 

-0.16543 ± j 1.49824* 

0.10975* 

3 

6.8825 

-0.29043 ± j 2.65129 

0.10S89* 

4 

8.3137 

-0.32120 ± j 2.91758 

0. 10943* 

5 

10.390 

-0.25878 ± j 3.22166 

0.08007 

6 

18.333 

-0.73126 ± j 4.20431 

0.17136 

7 

23.871 

-0.76931 ± j 4.79604 

0.15838 

S 

27.481 

-0.52110 ± j 5.18450 

0.10001 

9 

32.038 

-0.73261 ± j 5.56344 

0.13056 

10 

34.259 

-0.88063 ± j 5.78363 

0.15053 

11 

40.272 

-0.88313 ± j 6.29546 

0.13892 

12 

47.972 

-0.23151 ± j 6.78034 

0.03412 


* : critical constraints 



TABLE 18 


Natural Frequencies, Closed-Loop Eigenvalues and Modal Damping Ratios 
for Final Design of Case 4 in Example 2: ACOSS FOUR 


Mode 

Number 

Natural 

Frequency 

(w 2 ) 

Closed-loop 

Eigenvalue 

U i j w j) 

Damping 

Ratio 

ii) 

l 

1.7527 

-0.14811 ± j 1.33743* 

0.11007* 

"> 

2.1804 

-0.16354 ± j 1.49360* 

0.10884* 

3 

6.7730 

-0.29237 ± j 2.65644 

0.10940* 

4 

7.5292 

-0.30754 ± j 2.81009 

0.10879* 

5 

8.8137 

-0.41410 ± j 3.02121 

0.13579 

6 

16.516 

-0.74291 ± j 4.02792 

0.18138 

7 

23.311 

-0.77221 ± j 4.77106 

0.15977 

S 

24.868 

-0.69572 ± j 4.96686 

0.13872 

9 

27.885 

-0.65279 ± j 5.20579 

0.12442 

10 

30.156 

— 0.78313 ± j 5.34630 

0.14493 

1 1 

36.127 

-0.76335 ± j 5.S5341 

0.12932 

12 

47.773 

-0.24766 ± j 6.76245 

0.03660 


* : critical constraints 
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TABLE 19 


Iteration History for Example 2: ACOSS FOUR, Cases 1-4 

Total \Iass(*) 


Analysis 

Case 1 

Case 2 

Case 3 

Case 4 

1 

51.6970 

51.6970 

51.6970 

51.6970 

2 

45.2192 

47.0839 

50.0280 

25.6928 

3 

40.5850 

43.4828 

33.0527 

24.2318 

4 

29.9721 

32.6651 

31.6158 

23.3435 

5 

2.3.9877 

25.7153 

28.1675 

22.9081 

6 

22.8540 

22.4004 

26.2944 

22.9354 

7 

22.5623 

22.3012 

25.0935 

22.9354 

8 

22.4889 

22.2166 

24.1897 

22.9354 

4 

22.5179 

22.1845 

23.6746 


10 

22.5179 

22.1586 

23.6389 


11 

22.5179 

22.1404 

23.1769 


12 


22.1236 

23.1769 


13 



23.1769 



* : total mass Includes 8 units of nonstructural mass in addition to structural mass. 
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TABLE 20 


Iteration History for Example 3: Antenna Structure, Cases 1-10 
Control Design Variable Linking (I) 


Total Mass(kg) 


Analysis 

Case 1 

(144*) 

Case 2 
(72*) 

Case 3 
(36*) 

Case 4 
(18*) 

Case 5 
(8*) 

i 

502.14 

502.14 

502.14 

502.14 

502.14 

2 

462.20 

462.39 

466.52 

462.25 

455.96 

3 

345.34 

343.47 

349.71 

351.13 

359.45 

4 

254.17 

259.92 

261.19 

263.19 

283.10 

5 

213.91 

218.31 

235.33 

230.47 

231.57 

6 

184.11 

191.87 

210.26 

213.59 

212.83 

7 

174.34 

176.98 

199.19 

200.50 

204.21 

S 

168.93 

170.94 

192.61 

195.82 

200.82 

9 

169.35 

167.01 

193.88 

194.60 

198.91 

10 

165.61 

166.99 

190.74 

192.95 

198.71 

1 1 

164.42 

165.43 

188.11 

191.28 

198.45 


164.08 

165.03 

185.03 

189.94 

198.01 

m EM 

164.17 

164.84 

185.75 

189.03 

197.42 

MEM 

163.99 

164.64 

184.26 

187.51 

196.93 

'5 

163.76 

164.39 

182.27 

186.92 

196.46 


163.20 

164.19 

180.79 

186.66 

196.32 


163.10 

164.10 

179.90 

186.50 

196.31 


163.11 

164.05 

179.79 

186.34 


n 



179.64 




Analysis 

Case 6 
H*) 

Case 7 
(4*) 

Case 8 
(2*) 

Case 9 

a*) 

Case 10 

(i*) j 

1 

502.14 

502.14 

502.14 

502.14 

502.14 

2 

469.64 

465.14 

485.24 

474.45 

484.61 

3 

367.65 

361.18 

383.18 

374.55 

387.64 

4 

293.90 

280.72 

302.66 

291.97 

297.69 

5 

237.33 

231.87 

241.31 

234.26 

240.14 

6 

218.09 

212.22 

216.56 

214.15 

215.05 

7 

205.11 

204.82 

208.12 

207.01 

208.60 

8 

201.28 

201.05 

206.07 

207.00 

206.33 

9 

200.51 

200.69 

205.78 

204.60 

206.07 

10 

200.35 

200.63 

204.22 

204.47 

206.06 

11 

200.27 

200.59 

204.19 

204.46 

206.06 

12 



204.16 




* : number of independent control design variables 
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TABLE 21 


Final Cross Sectional Dimensions for Example 3: Antenna Structure, Cases 1-10 

Control Design Variable Linking (I) 


Pinal Cross Sectional Dimensions(cm) 


Case 

Plement 

1 

Element 

2 

Element 

3,4 

Element 

5,6 

Element 

7,8 


B 

25.00 b 

25.00 b 

23.24 

25.00 b 

19.31 

l 

H 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00 b 


T 

0.1000° 

0.1108 

0.1000° 

0.1935 

0.1000° 


B 

25.00 b 

25.00 b 

24.11 

25.00 b 

18.68 

o 

H 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

24.99 


T 

0.1046 

0.1129 

0.1000° 

0.1908 

0.1 OOF 


B 

25.00 b 

25.00 b 

17.68 

25.00 b 

18.87 

3 

H 

25.00 b 

25.00 b 

24.41 

25.00 b 

25.00 b 


1 

0.1393 

0.1263 

0.1000° 

0.2202 

0.1000° 


B 

25.00 b 

25.00 b 

16.67 

25.00 b 

18.56 

4 

H 

25.00 b 

25.00 b 

21.92 

25.00 b 

25.00 b 


T 

0. 1454 

0.1213 

0.1000° 

0.2490 

o.ioom 


B 

25.00 b 

25.00 b 

21.55 

25.00 b 

13.82 

S 

ii 

i 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00 b 


0.1748 

0.1152 

0.1000° 

0.2526 

0.1000° 


B 

25.00 b 

25.00 b 

22.04 

25.00 b 

10.84 

6 

II 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

24.83 


r 

0.2166 

0.1496 

0.1000° 

0.2134 

0.1000° 


B 

25.00 b 

25.00 h 

20.95 

25.00 b 

12.08 

7 

El 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

24.86 


I 

0.2054 

0.1360 

0.1000° 

0.2323 

0.1000° 


B 

25.00 b 

25.00 b 

19.04 

25.00 b 

14.43 

8 

H 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00 b 


T 

0.1790 

0.1115 

0.1000° 

0.2832 

0.1000° 


B 

25.00 b 

25.00 b 

18.47 

25.00 b 

14.50 

9 

II 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00° 


T 

0.1752 

0.1129 

0.1000° 

0.2883 

0.1000° 


B 

25.00 b 

25.00 b 

19.13 

25.00 b 

15.80 

10 

II 

T 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

0.1904 

0.1000“ 

0.1000° 

0.2815 

0.1000° 


u indicates lower bound value 


b indicates upper bound value 
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TABLE 22 


Summary of Example 3: Antenna Structure, Cases I - 1 0 
Control Design Variable Linking (I) 


Case 

Number( + ) 

Final Mass, ks> 

(**) 

Re{?.) 

Critical Constraints 

lmag(X) peak 

displacement 

peak 

control force 

1(144) 

163.11(18) 

a 6 ^10 


tflO ^7 1 6 

w i ^ «. 

- C2) 

164.05(18) 

- 

"4, 

<7io <7i« 

w, u 3 

3 (36) 

179.64(19) 



<7l0 #16 

W l W2 Wj ^ 

4(18) 

186.34(18) 

^4 6 10 


#10 #16 

^1 ^2 Wj ^ 

5 (8) 

196.31(17) 

<Tg <T 9 


#10 #16 

u \ w 3 K* 

6 (4) 

200.27(1 1) 

- 

l0 -u 

£ 7 i o f 7 1 6 

W| ^ 

7 (4) 

200.59(11) 

- 


#16 

w > w 3 K, 

S (2) 

204.16(12) 

- 

“V. 

#i 6 

w 3 K, 

9 (1) 

204.46(11) 

- 

W «4 

< 7 l 6 

w 3 

10 (1) 

206.05(10) 

- 

^4 

#10 #16 

^ 


+ : number of independent control design variables 

** : number of total analyses 
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TABLE 23 


Iteration History for Example 4: Antenna Structure, Linking on [//] 
Control Design Variable Linking (II) 


Total Mass(kg) 

Analysis 

Case 1 

Case 2 

Case 3 

Case 4 

Case 5 


(i*) 

(2*) 

(3*) 

(4») 

(5 + ) 

1 

502.14 

502.14 

502.14 

502.14 

502.14 

2 

484.61 

470.88 

471.87 

444.12 

454.83 

3 

387.64 

328.11 

328.47 

303.50 

346.62 

4 

297.69 

260.67 

254.76 

233.56 

268.66 

5 

240.14 

217.39 

214.99 

193.83 

215.70 

6 

215.05 

199.59 

200.57 

178.76 

185.24 

7 

208.60 

193.04 

192.60 

174.18 

175.46 

8 

206.33 

191.20 

190.47 

173.05 

176.56 

9 

206.07 

190.90 

190.23 

172.50 

173.25 

10 

206.06 

190.80 

190.14 

172.46 

172.66 

1 1 

206.06 

190.70 

189.94 

172.46 

172.46 

12 



189.85 


172.37 

13 



189.81 


172.32 



Analysis 

Case 6 

Case 7 

Case 8 

Case 9 

Case 10 


(6*) 

(7*) 

(8*) 

(9*) 

(10*) 

l 

502.14 

502.14 

502.14 

502.14 

502.14 

i 

453.14 

449.81 

450.11 

453.15 

469.10 

3 

299.54 

293.56 

294.62 

307.14 

350.45 

4 

228.72 

217.54 

218.41 

228.31 

268.09 

5 

187.80 

184.73 

180.21 

188.59 

215.51 

6 

175.65 

176.29 

174.64 

175.95 

181.84 


173.27 

173.07 

174.70 

173.41 

173.50 


172.41 

172.28 

174.38 

172.68 

171.19 


172.29 

171.93 

173.93 

171.10 

171.02 


172.17 

171.80 

171.27 

170.84 

170.96 



171.73 

170.97 

170.70 

170.92 




170.84 

170.65 


13 



170.81 




* : number of independent control design variables 
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TABLE 24 


Final Cross Sectional Dimensions for Example 4: 

Antenna Structure, Linking on [//], Control Design Variable Linking (II) 


Final Cross Sectional Dirncnsions(cm) 


c: 

ase 

Element 

1 

Element 

2 

Element 

3,4 

Element 

5.6 

Element 

7,8 


B 

25.00 b 

25.00 b 

19.62 

25.00 b 

15.51 

i 

H 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00” 


T 

0.1899 

0.1000° 

0.1000° 

0.2816 

0.1000° 


13 

25.00 b 

25.00 b 

19.78 

25.00 b 

24.94 

■"> 

H 

25.00 b 

25.00 b 

25.00 b 

23.05 

22.23 


T 

0.1830 

0.1326 

0.1000° 

0.2103 

0.1000° 


B 

25.00 b 

25.00 b 

20.03 

25.00 b 

22.86 

3 

H 

25.00 b 

25.00 b 

25.00 b 

23.37 

21.94 


T 

0.1869 

0.1413 

0.1000° 

0.2012 

0.1000° 


B 

25.00 b 

25.00 b 

20.90 

25.00 b 

20.51 

4 

H 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00” 


I 

0.1000° 

0.1000° 

0.1000“ 

0.2356 

0.1000° 


13 

25.00 b 

25.00*’ 

18.85 

25.00 b 

23.70 

5 

11 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00” 


T 

0.1025 

0.1036 

0.1000° 

0.2287 

0.1000° 


B 

25.00 b 

25.00 b 

18.72 

25.00 b 

21.03 

6 

II 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00” 


T 

0.1000° 

0.1107 

0.1000° 

0.2330 

0.1000° 


B 

25.00 b 

25.00 b 

17.65 

25.00 b 

20.47 

"7 

H 

25.00 b 

25.00 b 

24.89 

25.00 b 

25.00” 


T 

0.1000° 

0.1165 

0.1000° 

0.2318 

0.1000° 


B 

25.00 b 

25.00 b 

14.95 

25.00 b 

23.42 

s 

II 

25.00 b 

25.00 b 

23.00 

25.00” 

25.00” 


r 

0.1073 

0.1244 

0.1000° 

0.2203 

0.1000° 


B 

25.00 b 

25.00 b 

14.62 

25.00 b 

22.12 

9 

H 

25.00 b 

25.00 b 

22.36 

25.00 b 

25.00” 


T 

0.1013 

0.1319 

0.1000° 

0.2267 

0.1000° 


B 

25.00 b 

25.00 b 

19.96 

25.00 b 

20.29 

10 

H 

25.00 b 

25.00 b 

25.00 b 

25.00” 

25.00” 


T 

0.1000° 

0.1058 

0.1000° 

0.2296 

0.1000° 


Q indicates lower bound value 


b indicates upper bound value 
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TABLE 25 


Iteration History for Example 4: Antenna Structure, 
Linking on [// p ] and [//„], Control Design Variable Linking (II) 


Total Mass(kg) 


Analysis 

Case 11 

Case 12 

Case 13 

Case 14 

Case 15 


(2*) 

(4*) 

(6*) 

<»♦> 

(10*) 

1 

502.14 

502.14 

502.14 

502.14 

502.14 

2 

485.24 

482.76 

462.89 

420.35 

470.72 

3 

383.18 

343.98 

339.86 

288.81 

326.51 

4 

302.66 

276.53 

260.17 

229.95 

247.00 

5 

241.31 

225.43 

222.33 

196.72 

207.82 

6 

216.56 

200.56 

202.26 

180.59 

185.02 

7 

208.12 

189.44 

191.17 

174.55 

173.76 

8 

206.07 

187.41 

186.95 

172.76 

171.10 

9 

205.78 

186.31 

185.45 

171.64 

171.25 

10 

204.22 

1S5.54 

183.60 

170.30 

170.49 

11 

204.19 

185.21 

181.66 

170.41 

169.99 

12 

204.16 

185.13 

181.23 

170.37 

169.80 

13 


185.09 

180.70 


169.64 

14 



178.65 


169.55 

15 



178.06 



16 



177.86 



17 



177.31 



18 



176.72 



19 



176.67 






176.56 




♦ : number of independent control design variables 




TABLE 25 


Iteration History for Example 4: Antenna Structure. 

Linking on [// p ] and [// v ], Control Design Variable Linking (II), Continued 


Total Mass(kg) 


Analysis 

Case 16 
(12*) 

Case 17 
(14*) 

Case 18 
(16*) 

Case 19 
(18*) 

Case 20 
(20*) 

1 

502.14 

502.14 

502.14 

502.14 

502.14 

i 

475.80 

474.66 

474.47 

466.40 

434.64 

3 

301.89 

304.10 

304.66 

304.00 

287.18 

4 

228.23 

230.32 

231.61 

226.62 

224.61 

5 

188.16 

197.86 

197.41 

185.11 

183.32 

6 

175.46 

178.77 

178.20 

174.03 

172.69 

7 

174.21 

172.28 

172.25 

172.93 

171.31 

S 

173.71 

171.81 

171.09 

172.39 

169.42 

9 

172.95 

169.76 

170.49 

172.02 

166.81 

10 

171.19 

171.33 

170.21 

171.80 

166.79 

11 

170.91 

170.02 

169.72 

170.04 

166.22 

12 

170.69 

169.07 

169.16 

169.67 

165.66 

13 

169.68 

167.59 

167.78 

169.37 

165.61 

14 

169.49 

166.91 

167.26 

169.17 

166.50 

15 

168.38 

166.23 

167.00 

168.92 

165.30 

16 

167.40 

166.17 

166.65 

167.74 

165.14 

17 

167.01 

166.01 

166.29 

166.91 

164.95 

18 

166.85 


165.69 

165.84 

164.86 

19 

166.76 


165.46 

165.64 

164.80 

20 



165.34 

165.49 


21 



165.29 

165.38 



: number of independent control design variables 





TABLE 26 


Final Cross Sectional Dimensions for Example 4: 

Linking on [//,] and [//J, Control Design Variable Linking (II) 


Final Cross Sectional Dimensions(cm) 


Case 

Element 

1 

Element 

2 

Element 

3.4 

Element 

5.6 

Element 

7,8 

B 

11 II 
T 

25.00 b 

25.00 b 

0.1777 

25.00 b 

25.00 b 

0.1149 

18.99 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.2835 

14.13 

25.00 b 

0.1000° 

B 

12 II 
T 

25.00 b 

25.00 b 

0.1687 

25.00 b 

25.00 b 

0.1361 

19.57 

25.00 b 

0.1000° 

25.00 b 

22.84 

0.2032 

25.00 b 

22.14 

0.1000° 

B 

13 II 

T 

25.00 b 

25.00 b 

0.1239 

25.00 b 

25.00 b 

0.1259 

13.22 

25.00 b 

0.1000° 

25.00 b 

23.70 

0.2271 

25.00 b 

24.65 

0.1000° 

B 

14 H 
T 

25.00 b 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.1021 

22.11 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.2256 

20.00 

25.00 b 

0.1000° 

B 

15 II 
I 

25.00 b 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.1051 

19.45 

24.70 

0.1000° 

25.00 b 

25.00 b 

0.2246 

21.32 

25.00 b 

0.1000° 

B 

16 11 

T 

25.00 b 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.1341 

13.02 

24.12 

0.1000° 

25.00 b 

24.93 

0.2082 

24.32 

25.00 b 

0.1000° 

B 

17 II 
T 

25.00 b 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.1289 

13.64 

24.06 

0.1000° 

25.00 b 

25.00 b 

0.2061 

24.70 

25.00 b 

0.1000° 

B 

IS II 

r 

25.00 b 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.1213 

15.33 

23.14 

0.1000° 

25.00 b 

25.00 b 

0.2047 

25.00 b 

25.00 b 

0.1000° 

B 

19 II 
T 

25.00 b 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.1186 

14.99 

23.31 

0.1000° 

25.00 b 

25.00 b 

0.2084 

24.29 

25.00 b 

0.1000“ 

B 

20 H 

T 

25.00 b 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.1205 

14.91 

22.93 

0.1000° 

25.00 b 

25.00 b 

0.2054 

24.62 

25.00 b 

0.1000° 


u indicates lower bound value 


b indicates upper bound value 
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TABLE 27 


Summary of Example 4: Antenna Structure, Cases 1-20 
Control Design Variable Linking (II) 


Case 

Number(*) 

Pinal Mass, ka 
(**) 

Re{>.) 

Critical Constraints 

Imag(X) peak 

displacement 

peak 

control force 

i to 

206.06(11) 

- 


‘ 7.6 


2 ( 2 ) 

190.70(1 1) 

(T g <T 9 

a) i 4 & Js 

<7l6 


3(3) 

189.81(13) 

CT 4 (Tg (T 9 


<7l6 


4(4) 

172.46(1 1) 

G 3 0*4 ^8 


<7l6 


5 (5) 

172.32(13) 

^8 


<7l6 

Ut 

6(5) 

172.17(10) 



<7l6 


7(7) 

171.73(11) 


"■6 

<7l6 


S(S) 

170.81(13) 


<*>4, ^ 

<7l6 


9 (9) 

170.77(12) 

<J i <J S 

CD A CD j 

<7l6 

w 3 K* 

10 ( 10 ) 

170.92(11) 

- 

CD j 

*7 1 6 

w 3 ^ 

11 ( 2 ) 

204.16(12) 

- 


*6 

Wj ^ 

12(4) 

185.09(13) 

(Tg 

«., 5 

*7 16 

U* 

13(6) 

176.56(20) 

(J 4 <J g (T 9 


<7l6 

U* 

14 (S) 

170.37(12) 

^8 


<7l6 

u* 

15(10) 

169.55(14) 

(T 4 <T g 

"«4 

<7t6 

U* 

16(12) 

166.76(19) 

^3 ^4 ^8 


*16 

U 2 K, 

17(14) 

166.01(17) 



<7l6 

“3 ^ 

18(16) 

165.29(21) 


" 4 , "a 5 

<7l6 

Wi W 3 K, 

19(18) 

165.38(21) 

CS 4 (T 7 (J g 


<7l6 

w 3 K, 

20 ( 20 ) 

164.80(19) 

<r 7 

** 

3 

3 

?I6 

«, W 3 V* 


♦ : number of independent control design variables 

** : number of total analyses 
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TABLE 28 


Iteration History for Example 5: Antenna Structure, Additional Constraints 

Total Mass(kg) 


Analysis Case 1 Case 2 Case 3 


Case 1 

Control Cfforl 
Constraint 

Case 2 

Acceleration 

Constraints 

Case 3 

Control Effort 
& Acceleration 
Constraints 

502.14 

502.14 

502.14 

475.98 

446.92 

476.33 

369.44 

345.03 

375.87 

281.91 

305.57 

309.06 

242.58 

286.80 

286.28 

222.34 

277.24 

276.98 

216.91 

276.31 

273.76 

215.93 

276.01 

273.24 

214.42 

275.55 

272.95 

213.99 

275.23 

272.72 

213.71 

274.91 

271.76 


274.60 

270.65 


274.45 

270.50 


274.06 

270.44 


268.62 

270.38 


268.38 



268.03 



267.94 



267.88 






TABLE 29 


Final Cross Sectional Dimensions for Example 5: 
Antenna Structure, Additional Constraints 

Final Cross Sectional Dimensions(cm) 


Case 

Element 

I 

Flemcnt 

2 

Element 

3.4 

Flemcnt 

5,6 

Element 

7,8 

B 

25.00 h 

25.00 b 

16.94 

25.00 b 

18 55 

1 II 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

r 

0.2476 

0.1027 

0.1000° 

0.2452 

0.1000° 

B 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

2 II 

22.50 

22.56 

25.00 b 

23.92 

22.40 

T 

0.2935 

0.2656 

0.2432 

0.1255 

0.1552 

B 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

25.00 b 

3 H 

25.00 b 

22.15 

25.00 b 

24.11 

22.41 

T 

0.2841 

0.2793 

0.2405 

0.1278 

0.1532 


“ indicates lower bound value 


b indicates upper bound value 
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TABLE 30 


Iteration History for Example 5: Antenna Structure 
Re-solving 2x2 Riccati Equations and Truncation ot Gain Matrix 


Total Mass(kg) 


Analysis 

Case 4 

Re-solving 2x2 
Riccati Equations 

Case 5 
Truncation 
of Gain Matrix 

1 

502.14 

502.14 

? 

469.10 

469.10 

3 

366.71 

350.41 

4 

281.45 

268.57 

5 

220.83 

214.43 

6 

190.46 

181.03 

7 

177.59 

172.36 

8 

174.16 

170.64 

9 

174.02 

171.37 

10 

175.31 

171.25 

1 1 

175.05 

170.89 

12 

174.96 

170.85 

13 

174.79 

170.72 

14 

174.77 

170.71 

15 

174.74 

170.69 
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TABLE 31 


Final Cross Sectional Dimensions for Example 5: 
Re-solving 2x2 Riccati Equations and Truncation of Gain Matrix 


l inal Cross Sectional Dimensions(cm) 


Case 

Llcment 

1 

Clement 

2 

Element 

3,4 

Clement 

5,6 

Clement 

7,8 

B 

4 H 

r 

25.00 b 

25.00 b 

0.1025 

25.00 b 

25.00 b 

0.1000° 

16.58 

24.63 

0.1000° 

25.00 b 

25.00 b 

0.2446 

23.67 

25.00 b 

0.1000° 

B 

5 11 

T 

25.00 b 

25.00 b 

0.1000° 

25.00 b 

25.00 b 

0.1000° 

18.71 

24.33 

0.1000° 

25.00 b 

25.00 b 

0.2296 

23.11 

25.00 b 

0.1000° 


a indicates lower bound value 


b indicates upper bound value 



TABLE 32 


Final Closed-Loop Eigenvalues for Example 5: 

Re-solving 2x2 Riccati Equations and Truncation of Gain Matrix 

A , = <j, + jco dl (?„ %) 


Mode 

Number 

Case 10 
of Cxample 4 

Case 4 

Re-solving 2x2 
Riccati equations 

Case 5 
Truncation 
of Gain Matrix 

1 

-2.73 ± j 1.61 

-2.52 ± j 2.24 

-2.77 ± j 1.56 

2 

— 1.86 ± j 4.97 

— 1.54 ± j 5.01 

-1.20 ± j 5.04 

3 

-.683 ± j 23.7 

-.504 ± j 23.7 

-.739 ± j 23.6 

4 

-.519 ± j 50.1 

-.499 ± j 50.1 

-.528 ± j 50.1 

5 

-.866 ± j 64.7 

-.505 ± j 63.1 

-.906 ± j 63.2 

6 

-1.08 ± j 103. 

-.846 ± j 102. 

-1.34 ± j 102. 

7 

-.933 ± j 184. 

-.647 ± j 184. 

-.507 ± j 184. 

S 

-.623 ± j 241. 

-.501 ± j 238. 

-.514 + j 238. 

9 

-2.50 ± j 320. 

-.938 ± j 319. 

-2.98 ± j 320. 

10 

-2.25 ± j 324. 

-1.40 + j 325. 

-1.72 ± j 324. 

11 

.977 ± j 428. 

-.353 x 10" 2 ± j 428. 

.110 x 10' 3 ± j 427. 

12 

.076 + j 619. 

.127 x 10~ 2 ± j 617. 

-.462 x 10“ 6 ± j 618. 

13 

. 1S3 ± j 637. 

-.227 x 10~ 2 ± j 627. 

-.387 x 10- J ± j 626. 

14 

.203 ± j 681. 

-.425 x 10- 2 ± j 684. 

.395 x 10-* + j 681. 

15 

.121 ± j 1015. 

-.313 x 10 1 ± j 995. 

.438 x 10-’ ± j 1000. 

16 

-.121 ± j 1128. 

.131 x 10 2 ± j 1124. 

-.670 x 10* 6 ± j 1121. 

17 

.504 ± j 1268. 

.245 x I0~ J ± j 1294. 

.287 x lO' 6 ± j 1286. 

18 

.339 ± j 1375. 

-.378 x 10~ 2 ± j 1380. 

.971 x 10-" ± j 1378. 
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TABLE 33 


Iteration History for Example 5: Antenna Structure 
Variable Mass Design Elements 


Total Mass(kg) 



9 





TABLE 34 


Final Cross Sectional Dimensions for Example 5: Antenna Structure 
Variable Mass Design Elements 


Final Cross Sectional Dimensions(cm) 


Case 

Element 

1 

Element 

2 

Element 

3,4 

Element 

5,6 

Element 

7,8 

B 

6 H 

r 

25.00 b 

25.00 b 

0.1000“ 

25.00 b 

25.00 b 

0.1000“ 

18.26 

25.00 b 

0.1000“ 

25.00 b 

25.00 b 

0.2330 

22.80 

25.00 b 

0.1000“ 

B 

7 H 

r 

25.00 b 

25.00 b 

0.1000“ 

25.00 b 

25.00 b 

0.1580 

25.00 b 

25.00 b 

0.1000“ 

25.00 b 

25.00 b 

0.1728 

25.00 b 

25.00 b 

0.1000“ 

B 

8 H 
T 

25.00 b 

25.00 b 

0.1191 

25.00 b 

25.00 b 

0.1000“ 

13.98 

22.28 

0.1000“ 

25.00 b 

25.00 b 

0.2382 

21.38 

25.00 b 

0.1000“ 


« indicates lower bound value b indicates upper bound value 


TABLE 35 

Final Actuator Masses for Example 5: Antenna Structure 


Variable Mass Design Elements 
Actuator Mass(kg) 

Actuator 

C^ase 6 

Case 7 

Case 8 

Number 

c, = 4 kg l 8.5 N 

c, = 8 kg j 8.5 N 

c, = 2 kgj 8.5 N 

I 

2.700 

4.641 

1.489 

2 and 4 

4.006 

7.958 

2.001 

3 

3.920 

7.773 

1.974 
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TABLE 36 


Summary of Example 5: Antenna Structure, Cases 1-8 
Additional Problems 






Critical Constraints 


Case 

Number 

Final Mass, kg 
(**) 

Re{>.) 

I mag (a) 

peak 

displacement 

peak 

control force 

additional 

Case 10 
of 

Hxample 4 

170.92(11) 

’ 


<7l« 

w 3 K* 

NA 

1 

213.71(11) 

<7 1 (T s 


*7)0 <7l6 

- 

CE 

i 

267.88(19) 

<J 9 


<7l6 

^ w 3 i<, 

Or, 

3 

270.38(15) 


OJ di 

^7 16 

- 

CE On 

4 

174.74(15) 

G 3 ^4 G $ (7 ^ 


</l 6 


NA 

S 

170 69(15) 



*7l6 

"3 *<• 

NA 

6 

170.17(12) 

- 


‘7l6 

«<S 

™.4 

n 

/ 

180.37(12) 

O 3 ^4 G G % 

«N 4 

^7 16 

«3 


s 

165.03(12) 


UJ-, 

d <X 

<7l6 

"3 * 

™4 


** : number of total analyses 

NA : not applicable 
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TABLE 37 


Initial Complex Eigenvalues, Grillage Structure 




Mode 

Number 

Open- Loop 

Closed- Loop 

i 

-.0001 ± j 2.29 (.006) 

-.874 + j 2.31 (35.4) 

2 

-.001 ± j 7.05 (.018) 

-.431 ± j 7.06 (6.09) 

3 

-.007 + j 16.4 (.041) 

-1.87 ± j 17.3 (10.8) 

4 

-.009 ± j 18.6 (.047) 

-.762 ± j 17.6 (4.33) 

5 

-.018 ± j 27.1 (.069) 

— 1.36 ± j 27.2 (5.00) 

6 

-.039 ± j 39.5 (.099) 

-.878 ± j 39.9 (2.20) 

7 

-.040 ± j 39.9 (.100) 

-2.00 ± j 40.1 (4.97) 

8 

-.062 + j 49.7 (.124) 

-3.76 ± j 49.2 (7.64) 

9 

-.091 + j 60.2 (.150) 

-1.90 ± j 60.1 (3.15) 

10 

-.119 ± j 68.9 (.172) 

-2.01 ± j 69.0 (2.92) 

11 

-.135 + j 73.6 (.184) 

-2.58 ± j 74.0 (3.49) 

12 

-.146 ± j 76.5 (.191) 

-.840 ± j 75.6 (1.11) 

13 

-.231 ± j 96.0 (.240) 

-2.03 ± j 95.8 (2.11) 

14 

-.267 ± j 103. (.259) 

-1.39 ± j 103. (1.35) 

15 

-.281 ± ] 106. (.265) 

-2.26 ± j 106. (2.13) 

16 

-.316+ j 112. (.281) 

-1.41 ± j 112. (1.26) 

17 

-.353 ± j 119. (.297) 

-2.04 ± j 119. (1.72) 

18 

-.413 ± j 129. (.321) 

-1.66 ± j 128. (1.29) 

19 

-.505 ± j 142. (.355) 

-2.81 ± j 142. (1.98) 

20 

-.524 ± j 145. (.362) 

— 1.65 ± j 145. (1.14) 
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TABLE 38 


[independent Control Design Variables, Grillage Structure 



Case 1 

Case 2 

Case 3 

Case 4 

Case 5 

Case 6 

CDV's 

a, = *** = a* 

*i. 

a,,a 2 , 

*i.- 

Ot j , 

*i, .» ,a 19t 



* 2 = - = *25 

= ••• = Otjs 

a s = - = 

*10 = = *25 

^20 = = a 25 

Number 







of Indep. 

i 

2 

3 

5 

10 

20 

CDV's 








TABLE 39 

Final Structural Variables, Grillage Structure 


Depths (cm) 


Design 

Variable 

Number 

Case l 

Case 2 

Case 3 

Case 4 

Case 5 

Case 6 

1 

0.2849 

0.2768 

0.2510 

0.2392 

0.2557 

0.2466 

2 

0.1698 

0.1806 

0.1936 

0.1985 

0.1191 

0.1125 

3 

0.1000° 

0.1000° 

0.1000° 

0.1000“ 

0.1000° 

0.1000° 

4 

0.1953 

0.1973 

0.1769 

0.1788 

0.1908 

0.1560 

5 

0.1502 

0.1479 

0.1429 

0.1373 

0.1000° 

0.1000° 

6 

0.3656 

0.3626 

0.3670 

0.3691 

0.3701 

0.4128 

7 

0.2299 

0.2180 

0.2221 

0.2222 

0.1622 

0.1137 

8 

0.4724 

0.5021 

0.4951 

0.4719 

0.5880 

0.5573 


0 indicates lower bound value 
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TABLE 40 


Iteration Histories, Grillage Structure 


Total Mass (lb* see 2 1 In) 


Analysis 

Case l 

Case 2 

Case 3 

Case 4 

Case 5 

Case 6 


(1*) 

(2*) 

(3*) 

(5*) 

(10*) 

(20*) 


0.1294 

0.1294 

0.1294 

0.1294 

0.1294 

0.1294 


0.1448 

0.1433 

0.1440 

0.1442 

0.1442 

0.1468 


0.1332 

0.1327 

0.1369 

0.1353 

0.1352 

0.1366 

4 

0.1272 

0.1295 

0.1317 

0.1258 

0.1261 

0.1298 

5 

0.1231 

0.1256 

0.1238 

0.1222 

0.1257 

0.1223 

6 

0.1211 

0.1247 

0.1206 

0.1212 

0.1215 

0.1168 

7 

0.1211 

0.1200 

0.1185 

0.1197 

0.1182 

0.1163 

8 

0.1198 

0.1189 

0.1172 

0.1178 

0.1170 

0.1139 

9 

0.1192 

0.1189 

0.1166 

0.1170 

0.1145 

0.1114 

10 

0.1192 

0.1188 

0.1172 

0.1174 

0.1136 

0.1103 

11 

0.1194 


0.1173 

0.1173 

0.1126 

0.1103 


0.1191 


0.1172 

0.1171 

0.1111 

0.1089 

■EM 

0.1191 


0.1177 

0.1162 

0.1107 

0.1088 




0.1173 

0.1162 

0.1111 

0.1080 




0.1168 

0.1162 

0.1099 

0.1065 

HI 



0.1170 

0.1159 

0.1098 

0.1053 

17 



0.1169 

0.1158 

0.1098 

0.1047 

IS 



0.1169 

0.1159 

0.1093 

0.1044 

19 




0.1160 

0.1091 

0.1039 





0.1158 

0.1091 

0.1036 

21 




0.1157 

0.1093 

0.1042 

yy 




0.1157 

0.1083 

0.1041 

23 





0.1080 

0.1040 

24 





0.1082 

0.1039 

25 





0.1083 


26 





0.1082 


27 





0.1083 



* : number of independent control design variables 
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TABLE 41 


Final Closed-Loop Eigenvalues, Grillage Structure 


MAT, 0/ °) 


Mode 

Number 

Case 1 

Case 2 

Case 3 

I 

-1.58 ± j 4.85 (30.9) 

-1.62 + j 5.03 (30.7) 

— 1.89 ± j 5.14(34.4) 

2 

-.438 ± j 9.76 (4.49) 

-.321 ± j 9.93 (3.23) 

— .103 ± j 9.90 (1.04) 

3 

-2.26 + j 20.4 ( 1 1.0) 

-1.75 ± j 20.4 (8.53) 

-2.36 ± j 20.0(11.7) 

4 

-1.53 ± j 24.7 (6.17) 

-1.08 ± j 25.4 (4.25) 

-2.10 ± j 24.7 (8.50) 

5 

-2.53 ± j 30.5 (8.26) 

-2.08 ± j 31.0 (6.70) 

-2.58 ± j 29.9 (8.58) 

6 

-.415 ± j 38.9 (1.06) 

-.404 ± j 39.6 (1.02+) 

-.423 ± j 38.7 (1.09) 

7 

-2.24 + j 40.1 (5.59) 

-1.96 + j 39.9 (4.91) 

-2.62 ± j 39.2 (6.66) 

S 

-1.36 + j 50.7 (2.68) 

-1.87 ± j 51.5 (3.63) 

— 15.5 ± j 47.9 (30.7) 

4 

-15.4 + j 51.7 (28.6) 

-10.9 ± j 53.2 (20.1) 

-2.44 + j 51.2(4.75) 

10 

-2.49 ± j 55.1 (4.51) 

-1.98 ± j 54.9 (3.60) 

-3.14 ± j 52.3 (5.99) 

1 1 

-1.14 ± j 62.4 (1.83) 

-.674 ± j 63.6 (1.06) 

-.995 ± j 62.3 (1.60) 

n 

- .634 ± j 63.3 (1.00*) 

-.638 ± j 63.9 (1.00+) 

-.625 + j 62.3 (1.00+) 


-3.38 ± j 74.1 (4 55) 

-2.56 + j 75.0 (3.40) 

-2.54 + j 74.0 (3.42) 


-.944 ± j 93.1 (1.01*) 

-.934 + j 91.7 ( 1.02*) 

-.975 + j 88.5 (1.10) 

15 

-17.5 + j 97.6 (17.7) 

-11.8 ± j 100. (11.7) 

-20.7 ± j 95.4 (21.2) 

16 

-2.26 ± j 104. (2.16) 

-2.88 ± j 101. (2.85) 

-3.96 ± j 99.3 (3.99) 

17 

-3.28 ± j 105. (3.12) 

—4.61 + j 105. (4.40) 

-2.13 ± j 99.9 (2.13) 

18 

-2.67 ± j 106. (2.51) 

-1.08 ± j 105. (1.03+) 

-1.04 ± j 101. (1.03+) 

19 

— 1.21 ± j 120. ( 1 .01 + ) 

-1.24 ± j 121. (1.03+) 

— 1.20 ± j 119. (1.01+) 

:o 

-1.43 ±j 142. (1.01*) 

-1.49 ± j 140. (1.06) 

— 1.36 + j 135. (1.00+) 


* indicate critical damping ratio constraints (0.999 < < 1.03) 
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TABLE 41 


Final Closed-Loop Eigenvalues, Grillage Structure, Continued 


'-,= <+ +jOJj, (L. °°) 


Mode 

Number 

Case 4 

Case 5 

Case 6 

1 

-2.08 ± j 5.19 (37.1) 

-2.22 + j 5.92 (35.1) 

-1.85 ± j 6.56 (27.1) 

2 

-.100 + j 9.76 (1.02*) 

-.115 ± j 10.5 (1.09) 

-.112 + j 10.7(1.05) 

3 

-10.4 + j 19.9 (46.2) 

-1.16 ± j 22.2 (5.21) 

-.209 + j 20.3 (1.03*) 

4 

— 1.44 ± j 22.0 (6.55) 

—4.02 ± j 22.6 (17.5) 

-6.55 ± j 22.2 (28.3) 

5 

-2.71 ± j 30.0 (9.01) 

-1.92 ± j 28.6(6.70) 

-.724 + j 26.2 (2.77) 

6 

-.384 ± j 38.1 (1.01*) 

-.397 ± j 37.0 (1.07) 

-11.4 ± j 35.1 (30.9) 

7 

-2.01 ± j 38.8 (5.19) 

— 1.27 ± j 38.6(3.30) 

-.605 ± j 35.4 (1.71) 

8 

-13.7 ± j 48.3 (27.2) 

-12.0 ± j 45.8 (25.4) 

-.397 ± j 37.3 (1.07) 

9 

-2.24 ± j 50.5 (4.42) 

-.534 + j 48.3 (1.10) 

-.666 ± j 44.3 (1.50) 

10 

-3.07 ± j 51.7 (5.92) 

-.615 ± j 48.7 (1.26) 

-.479 ± j 47.9 (1.00*) 

i 1 

- .627 + j 61.8 (1 02*) 

-.525 ± j 52.2 (1.00*) 

-.496 ± j 48.4 (1.02+) 

12 

-.948 ± j 62.7 (1.51) 

-.648 ± j 56.1 (1.15) 

-.539 ± j 52.8 (1.02+) 

13 

-2.23 ± j 72.7 (3.06) 

-3.55 ± j 65.0 (5.45) 

— 1.93 ± j 54.6 (3.54) 

14 

-.881 ± j 86.6 ( 1 .02*) 

-.749 + j 72.7 (1.03*) 

-3.35 + j 66.1 (5.06) 

15 

— 19.4 ± j 96.7(19.7) 

-2.08 + j 77.5 (2.68) 

-37.0 ± j 66.3 (48.8) 

16 

-2.12 ± j 97.2 (2.19) 

-2.63 ± j 79.0 (3.32) 

-.742 ± j 71.9 (1.03+) 

17 

-3.81 ± j 98.4 (3.87) 

-12.6 + j 85.6(14.5) 

-7.59+ j 75.7 (9.97) 

18 

-.986 ± j 98.7 ( 1 00*) 

-1.62 + j 99.0 (1.63) 

-1.04 + j 90.5 (1.15) 

19 

-1.19 + j 117. (1.01+) 

-2.36 + j 110. (2.15) 

-1.35 ± j 95.1 (1.42) 


-1.40 ± j 135. (1.04) 

-1.33 ± j 133. (1.00*) 

-1.34 ± j 131. (1.02+) 


* indicate critical damping ratio constraints (0.999 < < 1.03) 
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Figure 2: Example 1 - Cantilever Beam 
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Final Dimensions (cm) 


o.io- 


0 . 00 - 


121 


Is 


1 ( 40 ) 2 ( 20 ) 3 ( 2 ) 4 ( 1 ) 5 ( 1 ) 

Case Number (No. of CDV’s) 


Figure 4: Final Cross Sectional Dimensions for Example 1 - Cantilever Beam, 

Cases 1-5, Control Design Variable Linking (1) 


K 






2 3 4 5 6 7 8 9 10 11 12 

Number of Analysis 


Figure 5: Iteration History for Example 1 - Cantilever Beam, 

Cases 6-9, Different Initial Feedback Gains 
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1 2 3 4 5 6 7 B 9 10 11 12 13 

Number of Analysis 


Figure 8: 


Iteration Histories for Example 2 - ACOSS FOUR, Cases 1 and 2 
(to j, > 1.341, (Oj2 > 1 -6, > 0.15 ) 




I I I 1 1 1 

8 9 10 11 12 13 

Analysis 


pie 2 - ACOSS FOUR, Cases 3 and 4 
1.5, 0.1093 < £,.< 0.1 1, i= 1,... ,4) 


a 






Initial Closed-Loop Eigenvalues 



Real Part 


Figure 12: Comparison of Initial Closed-Loop Eigenvalues for Example 3 - 

Antenna Structure 
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Imaginary Part 



Objective Mass (kg) 


220 



0 2 4 6 8 10 12 14 16 18 20 


Number of Analysis 


Figure 13: Iteration History for Example 3 - Antenna Structure, 

Cases 1-10, Control Design Variable Linking (I) 


9 











Final Dimensions 

B (cm) 

29 
20 
1 9 

1 O 


1 23 + 397 99 ^0 



H (cm) 



Element 


T (cm) 



Case Number (No. of CDV’s) 


Figure 14: Final Structural Dimensions for Example 3 - Antenna Structure, 

Cases 1-10, Control Design Variable Linking (1) 


150 


mum 




Final Mass (kg) 



Figure 15: Number of Independent CDV vs. Final Mass - Example 3, 

Antenna Structure, Cases 1-10, Control Design Variable Linking (I) 
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Objective Mass (kg) 


( )* : Number of independent Control Design Variables 



2 4 6 8 10 12 14 16 18 20 


Number of Analysis 


Figure 16: Iteration History for Example 4 - Antenna Structure, 

Cases I- 10, Control Design Variable Linking (II) 
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Final Dimensions 


B (cm) 

23 
20 
1 3 
1 O 
3 
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1 2 3 4 3 6 7 B O lO 

H (cm) 

29 
20 
1 3 
1 O 
3 - 
O- 

1 2 3 4- 3 9 7 S BIO 

T (cm) 

0.4-0 
0.30 
0.20 
0.1 O 
0.00 

1 2 3 4- 3 S 7 S BIO 

Case Number 

Figure 17: Final Structural Dimensions for Example 4 - Antenna Structure, 

Cases 1-10, Control Design Variable Linking (II) 
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( )* : Number of Independent Control Design Variables 



Figure 18: Iteration History for Example 4 - Antenna Structure, 

Cases 10-20, Control Design Variable Linking (II) 
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Final Dimensions 


B (cm) 



H (cm) 



Element 


T (cm) 



Case Number 


Figure 19: Final Structural Dimensions for Example 4 - Antenna Structure, 

Cases 10-20, Control Design Variable Linking (II) 


155 


□ □□□□ 



Final Mass (kg) 
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Figure 20: Number of Independent CDV vs. Final Mass - Example 4, 

Antenna Structure, Control Design Variable Linking (II) 
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Objective Mass (k 


Figure 



21: Iteration History for Example 5 - Antenna Structure, 

Additional Constraints 




Objective Mass (kg) 



4 5 6 7 8 9 10 II 12 13 14 15 16 

Number of Analysis 


Figure 22: Iteration History for Example 5 - Antenna Structure, 

Updating and Truncation of Gain Matrices 
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4 5 6 7 8 9 10 11 12 


Number of Analysis 


Figure 23: Iteration History for Example 5 - Antenna Structure, 

Variable Mass Design Elements 





Element 


Case Number 


Figure 24: Final Structural Dimensions for Example 5 - Antenna Structure, 

Additional Constraints and/or Special Features 
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GI5DH0 



# : Node Number 

* : Independent DV Number 
$ : Actuator Number 



Grillage Structure External Load 


Figure 25: Example 6 - Grillage Structure 


161 


Objective Mass (lb-sec /in) 



Number of Analysis 


Figure 26: Iteration History for Example 6 - Grillage Structure, 

Cases 1-6, Control Design Variable Linking (II) 









Depth (in) 



Depth (in) 



Design Variable Number 


Figure 27. Final Structural Dimensions for Example 6 - Grillage Structure, 
Cases 1-6, Control Design Variable Linking (II) 
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Final Mass (lb-sec 2 /in) 



Figure 28: Number of Independent CDV vs. Final Mass - Example 6, 

Grillage Structure, Cases 1-6, Control Design Variable Linking (II) 
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